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and  optimizing  TPV  cells  for  different  possible  spectra.  It  gives  results  for  GaSb  and  In- 
GaAs  cells  optimized  to  the  AMO  spectrum  which  closely  match  published  data  as  well 
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EXECUTIVE  SUMMARY 


A  major  limitation  in  space  exploration  is  the  availability  of  power.  While  satel¬ 
lites  stationed  around  Earth  can  produce  an  adequate  supply  of  power  using  solar  arrays, 
satellites  farther  from  the  Sun  require  unreasonably  massive  arrays  to  produce  very  little 
power.  In  order  to  accommodate  the  need  for  power  on  deep  space  missions,  it  is  neces¬ 
sary  for  satellites  to  carry  a  source  capable  of  providing  adequate  power  for  the  life  of  the 
mission.  This  is  currently  done  using  radioisotope  thermoelectric  generators  (RTGs). 
Current  RTGs  consist  of  a  heat  source,  Pu-238,  and  a  thermocouple  that  inefficiently 
converts  heat  into  electricity. 

In  order  to  improve  the  overall  efficiency  of  RTGs,  a  better  thermoelectric  con¬ 
verter  is  needed  to  replace  the  thermocouple.  This  thesis  examines  the  possible  use  of 
thermophotovoltaic  (TPV)  cells  in  RTGs.  Two  cells  were  developed  and  optimized  for 
the  spectrum  from  a  1300K  blackbody  which  simulates  the  spectrum  of  an  RTG  heat 
source. 

Current  TPV  cells  are  optimized  primarily  using  a  trial- and -error  method.  Cells 
are  built  and  tested  with  different  thicknesses  and  doping  levels  in  order  to  the  most  effi¬ 
cient  design.  This  method  is  very  expensive  and  time  consuming.  This  thesis  presents  a 
model  that  can  predict  the  output  of  a  cell  under  various  spectra.  The  model  can  easily  be 
changed  to  simulate  different  cell  thicknesses  and  doping  levels.  Additionally,  with  a 
comprehensive  set  of  material  parameters,  the  model  can  be  changed  to  simulate  cells  of 
different  materials. 

Hie  models  presented  in  this  thesis  were  built  using  the  Silvaco  Virtual  Wafer 
Fabrication  software  package.  This  software  package  has  been  shown  by  other  research¬ 
ers  as  a  capable  tool  for  modeling  advanced  solar  cell  designs,  but  there  is  no  record  of 
any  TPV  models  ever  being  designed. 

In  order  to  prove  the  capabilities  of  this  software  package,  an  initial  cell  was  de¬ 
signed  using  a  very  well  documented  material,  gallium  arsenide  (GaAs).  Voltage -current 
characteristics  and  frequency  response  data  were  recorded  from  this  model  and  compared 
with  experimental  data  from  a  similar  cell.  Once  the  model  was  verified,  more  exotic 
materials  could  be  examined.  Gallium  antimonide  (GaSb)  and  indium  gallium  arsenide 

xvii 


(InGaAs)  were  examined  for  TPV  cells  because  these  are  two  of  the  most  researched  ma¬ 
terials  and,  as  such,  there  is  a  significant  amount  of  experimental  data  to  compare  with. 

As  with  the  GaAs  cell  model,  the  GaSb  and  InGaAs  models  were  first  optimized 
for  the  AMO  spectrum  and  compared  to  experimental  results.  Again,  the  voltage- current 
relationship  and  frequency  response  of  the  cells  were  measured  and  found  to  be  very 
close  to  available  experimental  data.  Once  these  models  were  determined  to  be  accurate, 
they  could  be  optimized  for  a  different  spectrum. 

Programs  were  written  using  Matlab  to  run  iterations  of  Silvaco  model  changing 
layer  thicknesses  and  doping  levels  in  the  cells.  Other  programs  gathered  data  and  calcu¬ 
lated  and  recorded  open-circuit  voltage,  short  circuit  current,  maximum  power,  fill  factor, 
and  efficiency  of  each  cell  in  order  to  determine  the  most  efficient  design.  These  results 
were  also  compared  to  experimental  data  on  similar  cells  and  found  to  be  very  close. 


I.  INTRODUCTION 


As  space-based-  missions  take  both  manned  and  unmanned  vehicles  farther  from 
Earth,  the  need  for  a  better  source  of  power  becomes  more  important.  The  intensity  of 
light  from  the  sun  decreases  with  the  square  of  the  distance  from  it  making  solar  cells  an 
inadequate  source  of  power  on  distant  missions  because  the  area  of  cells  required  would 

be  tremendous.  As  an  example,  aim2  solar  panel  at  earth  would  have  to  be  25  m2  at 
Jupiter  and  almost  2000  m2  at  Pluto  to  yield  the  same  output  power.  Batteries,  another 
common  source  of  power  on  satellites  are  used  primarily  only  as  storage  devices  because 
they  are  only  capable  of  providing  a  few  hundred  watt  ■  hr/kg  .  This  means  that  the 
amount  of  batteries  required  to  provide  primary  power  for  a  long  mission  would  be  too 
massive  to  even  launch  into  space.  Current  deep  space  satellites  use  an  alternative  source 
of  power  known  as  a  radioisotope  thermoelectric  generator  or  (RTG). 

An  RTG  consists  of  a  heat  source  and  a  means  to  convert  that  heat  into  electricity. 
Existing  RTGs  accomplish  this  using  thermocouples.  Although  these  are  very  simple  and 
reliable  thermoelectric  converters,  they  operate  at  only  a  few  percent  efficiency.  In  order 
to  improve  overall  efficiency  of  an  RTG,  the  thermocouple  will  have  to  be  replaced  by  a 
more  capable  thermoelectric  converter. 

The  most  attractive  replacement  for  the  thermocouple  is  currently  the  thermo- 
photovoltaic  (TPV)  cell.  Research  has  been  ongoing  in  this  field  since  the  1950’ s,  but  the 
exotic  materials  necessary  for  high  efficiency  cells  has  only  been  recently  available. 

Here,  several  models  will  be  developed  and  optimized  for  both  the  sun’s  spectrum,  AMO, 
and  the  spectrum  from  a  1300K  blackbody  which  represents  well  the  spectrum  from  an 
RTG  heat  source. 

Due  to  the  fact  that  TPV  cells  work  on  the  same  principles  as  solar  cells,  these  are 
discussed  in  Chapter  n.  TPV  cells  are  discussed  in  Chapter  III.  In  Chapter  IV,  the  entire 
RTG  module  is  discussed  including  the  safety  precautions  necessary  when  using  a  radio¬ 
active  heat  source.  Chapter  V  discusses  the  software  packages  used  in  these  models  and 
Chapter  VI  gives  the  optimized  cell  characteristics  along  with  results,  and  Chapter  VII 
contains  conclusions  and  recommendations  for  future  work.  ATLAS  code  used  in  this 

thesis  is  in  Appendix  A  and  Matlab  code  is  in  Appendix  B. 
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II.  SOLAR  CELLS 


This  chapter  contains  the  basic  principles  of  solar  cells.  This  information  is  im¬ 
portant  to  this  thesis  because  thermophotovoltaic  cells  work  on  the  same  principles  as 
solar  cells. 


A.  HISTORY 


The  history  of  the  solar  cell,  also  known  as  a  photovoltaic  cell,  can  be  traced  all 
the  way  back  to  1839  when  the  then  nineteen-year-old  French  physicist  Antoine-Cesar 
Becquerel  first  discovered  the  photovoltaic  effect  [Ref.  1].  While  experimenting  with 
electrolytic  cells  containing  two  different  metal  electrodes  in  a  solution,  he  found  that 
some  combinations  of  metal  and  solution  would  produce  a  small  current  when  exposed  to 
light. 

The  first  actual  solar  cell  was  built  in  1877  by  Charles  Fritts  [Ref.  1].  He  applied 
a  very  thin  gold  film  to  a  piece  of  selenium.  Although  this  cell  operated  with  less  than 
one  percent  efficiency,  it  was  able  to  produce  a  measurable  amount  of  power  [Ref.  1].  A 
similar  cell  was  constructed  in  1927  consisting  of  a  copper  film  on  piece  of  copper  oxide 
[Ref.  1].  This  cell  also  operated  with  less  than  1%  efficiency  [Ref.  1]. 

The  next  breakthrough  came  in  1941  when  Russell  Ohl  built  the  first  silicon  solar 
cell  [Ref.  4].  In  1954,  Pearson,  Chapin,  and  Fuller  expanded  on  this  design  and  were  able 
to  build  a  silicon  cell  with  an  efficiency  of  6%  [Ref.  1].  Since  that  time,  there  has  been  a 
great  deal  of  progress  in  solar  cell  design.  Gallium  arsenide  cell  designs  have  been  able 
to  achieve  efficiencies  of  over  20%  and  multijunction  cells  have  been  built  with  up  to 
37%  efficiency  [Ref.  1]. 
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B. 


THE  P-N  JUNCTION 


N-type  materials  have  electrons  as  the  majority  carriers  while  p-type  materials 
have  holes  as  the  majority  carriers.  What  is  important  to  semiconductor  physics  is  what 
happens  when  these  two  types  of  materials  come  together  [Ref.  2].  As  can  be  seen  from 
Figure  1  below,  when  the  two  materials  come  in  contact,  some  electrons  move  across  the 
junction  to  fill  holes  on  the  p-side  and  some  holes  move  across  the  junction  to  the  n-side 
[Ref.  2], 
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Figure  1.  p-n  Junction 


When  these  electrons  and  holes  migrate  across  the  junction,  an  electric  field 
forms.  This  area  is  referred  to  as  the  space  charge  region  or  depletion  region,  and  it  cre¬ 
ates  a  built  in  potential  in  the  material.  The  energy  band  diagram  in  Figure  2  [Ref.  2] 
helps  to  give  a  better  idea  what  happens  when  n-type  and  p-type  materials  are  brought 
together. 
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Figure  2.  Energy  Band  Model  of  a  p-n  Junction  [From  Ref.  2.] 

The  Fermi  level  is  constant  throughout  the  material,  so  a  straight  line  is  drawn  and 
the  two  band  diagrams  are  shifted  to  match  it.  In  order  to  have  continuous  conduction 
and  valence  bands  in  the  material,  the  bands  of  the  p-  and  n-type  material  must  bend  to 
meet  each  other  [Ref.  2].  This  bend  in  energy  bands  creates  a  built-in  potential  in  the  ma¬ 
terial  equal  to  Vbi  as  seen  in  Figure  3  [Ref.  2].  It  is  this  potential  that  causes  the  voltage 
drop  across  a  diode  [Ref.  2]. 


Figure  3.  Built-In  Potential  [After  Ref.  2.] 
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Figure  4  shows  how  the  junction  reacts  when  a  bias  is  placed  across  it  [Ref.  2]. 
Part  (a)  shows  the  band  diagram  under  equilibrium  conditions  with  no  bias  present.  This 
diagram  was  previously  explained. 

Part  (b)  shows  the  junction  under  forward  bias  conditions.  This  is  when  the  anode 
is  held  at  a  higher  potential  than  the  cathode.  As  more  voltage  is  applied  across  the  junc¬ 
tion,  it  acts  against  the  built  in  potential  Vbi .  Once  the  bias  across  the  junction  reaches  a 
level  that  is  equal  to  Vbi ,  current  will  begin  to  flow  through  the  junction.  The  diode  will 
then  act  like  a  short  circuit  with  a  constant  voltage  drop  of  Vbi . 

Part  (c)  shows  the  band  diagram  of  the  p-n  junction  under  reverse  bias  conditions. 
This  occurs  when  the  voltage  at  the  cathode  is  higher  than  the  voltage  at  the  anode.  In 
this  case,  the  potential  across  the  junction  acts  to  increase  the  level  between  the  energy 
bands  of  the  two  materials  thereby  preventing  current  from  flowing. 


6 


- Ev 

(a)  Equilibrium  (Va— 0) 


Ec 
Efn 
Ei 

Ev 

(b)  Forward  Bias  (Va>0) 


(c)  Reverse  Bias  (Va<0) 

Figure  4.  p-n  Junction  Under  Bias  [After  Ref.  2.] 
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c. 


PHOTOVOLTAIC  CELL 


A  photovoltaic  cell  [Ref.  3],  shown  in  Figure  5,  produces  power  using  photogen¬ 
eration.  When  a  photon  of  light  with  energy  higher  than  the  bandgap  energy  strikes  the 
material,  an  electron- hole  pair  is  created.  The  cell  itself  is  simply  a  p-n  junction  which 
was  discussed  in  the  previous  section.  The  electric  field  in  the  depletion  region  at  the 
junction  attracts  the  free  electrons  towards  the  cathode,  the  top  layer  in  the  figure,  and 
attracts  the  holes  towards  the  anode,  the  bottom  layer. 


Ugfct  Can«'»u«  ugM  »  Absorbed 

Emetrci  end  Hoi*  at  Back 


O  C'kClron 
©  Holt 


Figure  5.  Simple  Solar  Cell  [After  Ref.  3.] 


Current  is  carried  away  from  the  cell  by  contacts  which  are  connected  both  on  the 
top  and  bottom  of  the  cell.  The  contact  on  the  top  of  the  cell  is  typically  made  up  of  a 
thin  metal  electrical  grid  to  avoid  shading  the  cell.  The  bottom  metal  contact,  on  the 
other  hand,  spans  the  entire  surface  of  the  cell.  These  contacts  are  connected  to  the  ex¬ 
ternal  load. 

D.  INCIDENT  SPECTRUM 

Solar  cells  can  be  designed  to  be  optimized  for  different  spectra  based  on  their  in¬ 


tended  use  [Ref.  3].  Cells  intended  for  outer  space  are  designed  to  the  spectrum  of  Air 
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Mass  Zero  (AMO).  This  is  the  spectrum  from  the  sun  measured  just  outside  the  earth’s 
atmosphere.  Figure  6  shows  the  AMO  spectrum  in  |^W/m2/pmJ.  The  data  was  meas¬ 
ured  by  The  National  Renewable  Energy  Laboratory  (NREL)  [Ref.  4].  The  total  power 
in  this  spectrum  is  approximately  1367  W/m2  . 


As  mentioned  previously,  the  energy  of  a  photon  required  to  cause  photogenera¬ 
tion  in  a  material  must  be  greater  than  the  bandgap  energy  of  the  material.  Because  of 
this  relationship,  it  is  often  useful  to  see  the  AMO  spectrum  as  a  function  of  energy  in¬ 
stead  of  wavelength.  The  relationship  of  photon  energy  to  wavelength  is  expressed  as 

he 


E  = 


X 


(2.1) 


E  [eV  ] 


1.24 

?i[|tm] 


(2.2) 
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where  h  is  Planck’s  constant,  c  is  the  speed  of  light,  and  X  is  the  wavelength  of  light 
[Ref.  3].  Figure  7  shows  the  AMO  spectrum  [Ref.  4]  as  a  function  of  energy. 


As  the  sun’s  spectrum  passes  through  the  earth’s  atmosphere  some  wavelengths 
are  absorbed  and  others  are  reflected.  This  causes  the  radiance  at  the  earth’s  surface  to  be 
a  subset  of  the  AMO  spectrum.  Other  common  spectra  are  Air  Mass  One  (AMI)  and  Air 
Mass  One  Point  Five  (AMI. 5).  AMI  is  the  spectrum  from  the  sun  at  earth’s  surface  on 
the  equator  when  the  sky  is  clear.  Moving  farther  from  the  equator  to  higher  latitudes 
and  the  presence  of  clouds  and  other  weather  conditions  further  reduces  the  radiance  of 
the  spectrum  that  reaches  earth’s  surface.  AMI. 5  is  commonly  used  for  testing  terrestrial 
based  solar  cells  [Ref.  3]. 
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E. 


POWER  OUTPUT 


As  discussed  previously,  a  photovoltaic  cell  is  simply  a  p-n  junction.  It  therefore 
has  associated  with  it  the  same  current  voltage  relationship  as  a  diode  when  there  is  no 
light  incident  on  the  cell. 

Figure  8  shows  currents  that  contribute  to  the  total  current  from  a  solar  cell  [Ref. 
5].  The  curve  in  part  (a)  shows  the  photocurrent  from  the  cell.  This  is  a  constant  current 
level  that  is  dependent  on  the  amount  power  incident  on  the  cell.  The  dark  current  curve 
in  part  (b)  is  an  exponential  curve  that  is  a  characteristic  of  the  material.  To  find  the  cur¬ 
rent- voltage  relationship  of  the  solar  cell,  the  dark  current  must  be  subtracted  from  the 
photocurrent.  Part  (c)  of  Figure  8  shows  the  resulting  I-V  curve  of  the  cell. 


(a)  Photocurrent 


Figure  8. 


(b)  Dark  Current  (c)  |.y  Curve 

Current  in  a  Solar  Cell  [After  Ref.  5] 


Figure  9  shows  an  example  of  a  solar  cell  I-V  curve.  The  points  on  this  curve  of 
most  importance  are  labeled.  When  the  anode  and  cathode  are  shorted,  the  maximum 
current,  called  short  circuit  current,  flows  through  the  circuit  [Ref.  3].  This  is  labeled  as 
Isc  on  the  figure.  If  the  circuit  is  opened,  no  current  will  flow;  however  the  maximum 
voltage,  called  open-circuit  voltage,  can  be  measured  across  the  cell  [Ref.  3]. 
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While  the  cell  can  operate  anywhere  on  this  curve,  it  is  most  efficient  when  oper¬ 
ated  at  the  maximum  power  point,  labeled  Pma  in  the  figure  [Ref.  3].  This  is  the  point  on 
the  curve  where  the  product  of  the  operating  voltage  and  current  is  the  highest.  The  effi¬ 
ciency  of  the  cell  is  calculated  as  [Ref.  3] 

n  __™^,  (2.3) 

p 

in 

where  Pm  is  the  power  in  the  incident  spectrum;  for  AMO  this  is  1367  W/m2  [Ref.  3]. 

Another  important  characteristic  of  a  solar  cell  is  the  fill  factor,  which  indicates  the 
squareness  of  the  I-V  curve  [Ref.  3].  This  is  calculated  as  [Ref.  3] 

V  I 

FF  =  mp  mp  .  (2.4) 

^ 0(1  sc 

The  dotted  box  that  intersects  I sc  and  Vrx:  in  the  figure  represents  the  optimal  I-V 
curve  with  a  fill  factor  of  one.  The  fill  factor  is  a  measure  of  how  much  of  this  dotted 
box  the  actual  FV  curve  occupies  [Ref.  3]. 
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F. 


EFFICIENCY  IMPROVEMENTS 


There  are  many  factors  that  can  be  changed  to  help  improve  the  efficiency  of  solar 

cells. 

1.  Reflection 

An  untreated  solar  cell  can  reflect  as  much  as  36%  of  the  incident  light  [Ref.  3]. 
Anti- reflective  coatings  such  as  silicon  nitride  can  reduce  the  reflectivity  of  the  cell  to  as 
low  5%  [Ref.  3].  In  addition  to  adding  antireflective  coatings,  texturing  the  cell  surface 
also  reduces  the  amount  of  light  reflected  at  the  surface  of  the  cell  [Ref.  3]. 

Figure  10  shows  an  example  of  surface  texturing  [Ref.  6].  This  cell  was  able  to 
operate  at  24%  efficiency  which  is  very  high  for  a  silicon  cell  [Ref.  6].  The  inverted 
pyramids  labeled  in  the  figure  are  at  such  an  angle  so  that  light  that  is  reflected  from  the 
surface  of  the  cell,  is  reflected  in  such  a  way  that  it  will  hit  one  of  the  other  pyramids 
where  it  may  be  absorbed  [Ref.  6].  The  high  absorption  of  this  cell  led  to  its  record  set¬ 
ting  efficiency. 


finger  "inverted"  pyramids 


rear  contact  oxide 


Figure  10.  Anti- Reflective  Coating  and  Textured  Cell  Surface  [From  Ref.  6.] 
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2.  Low  Energy  Photons 

Only  photons  with  energy  higher  than  the  bandgap  energy  of  the  material  can  cre¬ 
ate  an  electron- hole  pair  in  the  cell.  Photons  with  too  little  energy  to  create  an  electron- 
hole  pair  will  be  absorbed  and  increase  the  temperature  of  the  cell  [Ref.  3].  This  in  turn 
reduces  the  efficiency  of  the  cell.  Some  designs  utilize  optical  filters  to  block  low  energy 
photons  to  avoid  this  problem  [Ref.  3]. 

3.  High  Energy  Photons 

Photons  with  energy  significantly  higher  than  the  bandgap  energy  can  still  create 
an  electron- hole  pair,  but  the  excess  energy  is  given  off  as  heat  in  the  cell,  thus  reducing 
the  efficiency  [Ref.  3].  Some  of  these  photons  may  simply  pass  through  the  cell  without 
creating  an  electron- hole  pair.  A  back  surface  reflector  can  reflect  these  high  energy  pho¬ 
tons  back  through  the  cell  where  they  once  again  have  the  opportunity  to  create  an  elec¬ 
tron-hole  pair  [Ref.  3]. 


4.  Cell  Resistance 

There  are  several  sources  of  resistance  in  a  solar  cell.  There  is  a  resistance  in  the 
material  itself  that  impedes  the  path  of  free  electrons  and  holes  to  the  contacts  [Ref.  3]. 
There  are  also  junction  resistances  where  the  electrical  contacts  are  joined  to  the  cell 
[Ref.  3].  Finally  there  is  resistance  in  the  contacts  [Ref.  3].  All  of  these  resistances  cause 
losses  which  can  be  significant  especially  in  high  current  applications. 

5.  Material  Defects 

Defects  in  the  crystal  structure  aid  in  trap  assisted  recombination  which  reduces 
cell  efficiency  [Ref.  2].  This  is  a  significant  source  of  power  loss  in  inexpensive  poly¬ 
crystalline  and  amorphous  cells  reducing  efficiency  by  as  much  as  15%  [Ref.  3].  There 
are  also  a  significant  number  of  material  defects  when  lattice  mismatch  is  a  factor  such  as 
in  multijunction  cells.  The  effects  of  this  in  single  junction  cells  can  be  reduced  by  using 
high  quality  crystals  [Ref.  3].  The  fewer  the  defects  in  the  crystal  structure,  the  higher 
the  efficiency  of  the  cell.  In  multijunction  cells,  a  great  deal  of  research  has  been  done 
with  ternary  and  quaternary  compounds  (compounds  with  three  or  four  elements  respec- 
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tively)  [Ref.  3].  By  altering  the  molar  concentrations  of  the  elements  in  these  com¬ 
pounds,  the  lattice  constant  can  be  set  to  a  desired  value  reducing  the  amount  of  lattice 
mismatch. 


6.  Shading 

In  order  to  collect  free  carriers  from  the  cell,  there  must  be  metal  contacts  on  the 
surface.  These  contacts,  while  providing  a  necessary  function,  also  shade  the  parts  of  the 
cell  they  cover  from  incident  light  [Ref.  3].  Larger  and  more  closely  spaced  contacts 
have  the  benefit  of  being  able  to  collect  more  generated  carriers  and  having  lower  internal 
resistance  [Ref.  3].  Conversely  they  shade  more  of  the  cell  causing  fewer  carriers  to  be 
generated.  Smaller  and  more  loosely  spaced  carriers  shade  less  of  the  cell,  but  also  have 
higher  resistance  [Ref.  3].  In  addition  when  contacts  are  more  loosely  spaced  on  the  sur¬ 
face  of  the  cell,  carriers  are  more  likely  to  recombine  before  they  can  be  collected  by  the 
contacts  [Ref.  3].  By  trial  and  error,  it  has  been  found  that  the  optimal  contact  arrange¬ 
ment  consists  of  equally  spaced  contacts  covering  approximately  8%  of  the  cell’s  surface 
[Ref.  3], 

This  chapter  reviewed  the  basic  principles  of  solar  cells.  The  input  spectrum  and 
output  power  were  covered  along  with  methods  for  improving  the  efficiency.  The  next 
chapter  will  introduce  a  related  concept,  the  thermophotovoltaic  cell. 
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III.  THERMOPHOTO  VOLTAIC  CELLS 


This  chapter  discusses  the  basic  principles  of  thermophotovoltaic  (TPV)  cells. 
Thermophotovoltaic  cells  are  designed  to  convert  the  infrared  (IR)  region  of  the  electro¬ 
magnetic  spectrum  into  electricity.  The  goal  of  this  thesis  is  to  model  and  optimize  these 
cells,  so  it  is  important  to  know  some  background  information  on  them. 

A.  HISTORY 


Although  the  ‘Father’  of  thermophotovoltaics  is  unknown,  the  field  can  be  traced 
at  least  as  far  back  as  1956  to  Dr.  Henry  H.  Kolm  [Ref.  7].  While  working  at  MIT’s  Lin¬ 
coln  Laboratory,  Dr.  Kolm  built  a  primitive  TPV  device  using  a  Coleman  camping  lan¬ 
tern  as  a  radiator  and  silicon  solar  cells  for  photoconverters.  Although  his  system  only 
produced  one  watt  of  power,  it  was  a  start.  Throughout  the  1960’s,  researchers  at  MIT 
continued  research  on  TPV  devices  focusing  on  germanium  cells  [Ref.  7.].  In  the  1970’s 
thermophotovoltaics  was  studied  as  a  way  to  provide  a  portable  power  source  for  the  U.S. 
Army  [Ref.  7].  Early  TPV  generators  consisted  of  standard  silicon  or  germanium  solar 
cells  surrounding  a  heat  source  [Ref.  7].  Around  1980  a  general  interest  in  the  area  was 
lost  because  of  the  exceedingly  low  efficiency  of  these  early  TPV  cells.  With  the  recent 
advancements  in  III-V  semiconductor  technology  (materials  made  of  a  combination  of 
group  III  and  group  V  elements  from  the  periodic  table),  there  was  a  possibility  for  sig¬ 
nificantly  higher  efficiency  cells  [Ref.  7].  Recent  research  has  shown  that  gallium  anti- 
monide  (GaSb)  and  indium  gallium  arsenide  (InGaAs)  are  two  of  the  most  promising  ma¬ 
terials  for  new  TPV  cells. 
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B.  INCIDENT  SPECTRUM 


1.  Blackbody  Radiation 

Objects  with  a  temperature  abo\e  zero  Kelvin  radiate  energy.  In  order  to  estimate 
the  amount  of  energy  radiated,  it  is  convenient  to  introduce  the  concept  of  a  blackbody. 

A  blackbody  is  characterized  by  three  key  properties  [Ref.  8] : 

1.  100%  absorption  of  all  incident  radiation 

2.  Maximum  possible  emission  at  any  temperature  and  wavelength 

3.  Diffuse  emitter,  radiates  energy  equally  in  all  directions 


An  ideal  emitter,  a  blackbody,  radiates  energy  over  the  whole  electromagnetic 
spectrum  as  a  function  only  of  temperature  and  wavelength.  A  blackbody  emitter  radi¬ 
ates  energy  per  wavelength  into  a  unit  solid  angle  in  a  distribution  described  by  Planck’s 
Law  for  spectral  radiance  as  [Ref.  8] 
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where  h  is  Planck’s  constant,  c  is  the  speed  of  light  in  a  vacuum,  and  k  is  Boltzmann’s 
constant.  T  and  X  are  the  only  two  variables  and  refer  to  the  blackbody  temperature  in 
Kelvin  and  the  wavelength  of  energy  emitted  in  meters.  The  units  of  spectral  radiance 
are  bracketed. 


As  can  be  seen  from  Planck’s  Law  for  spectral  radiance,  at  a  given  temperature, 
the  amount  of  radiation  varies  as  function  of  wavelength.  An  important  feature  is  that  at 
any  wavelength,  the  magnitude  of  radiant  energy  increases  as  temperature  increases  [Ref. 
8].  This  means  that  a  hotter  blackbody  will  emit  more  energy  than  a  colder  one. 

Another  important  feature  of  Planck’s  Law  is  the  distribution  of  radiant  energy. 
More  radiant  energy  is  emitted  at  shorter  wavelengths  for  hotter  blackbodies  [Ref.  8]. 
This  phenomenon  is  described  by  Wien’s  Displacement  Law  [Ref.  8], 
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This  equation  gives  the  wavelength  at  which  the  maximum  radiation  is  emitted 
from  a  blackbody  at  a  given  temperature.  It  can  easily  be  seen  that  the  location  of  this 
peak  is  inversely  proportional  to  the  temperature  of  the  blackbody. 

In  order  to  find  the  total  power  emitted  by  a  blackbody,  it  is  necessary  to  integrate 
Planck’s  Law  over  all  wavelengths.  This  calculation  yields  the  Stefan-Boltzmann  Law 
[Ref.  8], 

Eh=aT4.  (3.3) 

where  o  =  5.670x10  ^W/m2  /K4J  .  Here  it  can  be  seen  that  the  total  radiated  power  is 

dependent  only  on  temperature. 


When  working  with  solar  cells,  the  spectrum  of  interest  is  AMO.  This  spectrum 
corresponds  to  the  radiation  from  blackbody  at  a  temperature  of  5800K.  Because  energy 
is  radiated  outward  from  a  body  at  an  average  angle  perpendicular  to  its  surface,  the  mag¬ 
nitude  of  radiation  in  a  given  direction  is  equal  to  only  one-fourth  the  total  radiation  [Ref. 
8].  This  corresponds  to  n  steradians.  With  this  information,  the  equation  for  the  hemi¬ 
spherical  spectral  radiant  flux  in  any  direction  from  a  flat  surface  emitter  can  be  derived 
as  [Ref.  8] 
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In  order  to  find  the  AMO  spectrum,  there  is  one  more  factor  that  must  be  consid¬ 
ered.  The  magnitude  of  radiation  emitted  from  a  source  is  inversely  proportional  to  the 
square  of  the  distance  from  the  source.  The  factor  used  to  calculate  the  AMO  spectrum  is 
found  by  taking  the  square  of  the  ratio  of  the  sun’s  radius  to  the  mean  distance  between 
the  earth  and  the  sun.  This  factor,  referred  to  as  the  dilution  factor,  comes  out  to 

2.165  xlO °  [Ref.  8].  The  final  equation  for  the  AMO  spectrum  can  now  be  expressed  as 
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where  /  is  the  dilution  factor. 


Using  this  equation,  Figure  11  below  shows  the  spectrum  from  a  5800K  black- 
body  adjusted  by  the  dilution  factor.  It  is  plotted  on  top  of  the  actual  AMO  spectrum  as 
measured  by  NREL  [Ref.  4] .  It  can  be  seen  that  the  two  spectra  are  very  closely  corre¬ 
lated  in  shape,  magnitude,  and  area. 


Figure  11.  AMO  vs.  5800K  Blackbody  Spectrum 


From  Planck’s  distribution  law,  the  magnitude  of  radiation  at  a  given  wavelength 
increases  with  temperature,  so  the  magnitude  of  radiation  from  the  sun  is  much  larger 
than  the  magnitude  of  radiation  from  any  heat  source  with  a  lower  temperature.  Unlike 
photovoltaic s,  however,  where  the  radiant  source,  the  sun,  is  millions  of  miles  away,  a 
heat  source  used  for  thermophotovoltaic  applications  can  be  placed  in  very  close  prox¬ 
imity  to  the  cell.  Although  the  sun  radiates  much  more  energy  than  a  heat  source,  be¬ 
cause  of  the  distance  between  the  source  and  the  cell,  the  intensity  of  radiation  at  the  cell 
is  much  higher  for  a  heat  source. 
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2.  Radioisotope  Thermoelectric  Generator  Spectrum 

When  dealing  with  solar  cells,  the  input  spectrum  is  a  subset  of  the  spectrum  from 


the  sun  [Ref.  3].  With  TPV  cells,  the  input  spectrum  is  dependent  on  the  temperature  of 
the  radiator.  Table  1  shows  several  temperatures  and  the  wavelength  at  which  their  radia¬ 
tion  spectrum  is  a  maximum  as  calculated  from  Equation  3.2. 


T(K)  | 

?Pe<JVm) 

300  i 

9.66 

1000  i 

2.90 

1300 

_ J_ 

2.23 

2000 

1.45 

3000 

0.97 

5000  ! 

i 

0.58 

5800  | 

0.50 

7000  j 

0.41 

10000  j 

0.29 

Table  1.  Temperature  vs.  Peak  Wavelength 


The  previous  section  discussed  blackbody  radiation  and  compared  a  5800K  black- 
body  radiation  curve  to  the  AMO  spectrum  in  Figure  12.  Because  these  two  were  so 
closely  correlated,  it  was  assumed  that  a  blackbody  curve  would  be  a  good  best  case  ap¬ 
proximation  for  the  input  spectrum  from  a  radioisotope  thermoelectric  generator  (RTG) 
heat  source. 

Current  RTG’s  use  Pu-238  as  a  heat  source.  In  a  pure  form,  Pu-238  can  reach 
surface  temperatures  of  1300K  [Ref.  9].  For  simulations  in  this  research,  the  input  spec¬ 
trum  was  assumed  to  be  a  1300K  blackbody.  Figure  24  shows  the  spectrum  that  was 
used  in  simulations. 
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Figure  12.  1300K  Blackbody  Spectrum 


C.  THERMOPHOTO V OLT AI C  SYSTEM 


There  are  more  parts  to  a  thermophotovoltaic  system  than  just  the  TPV  cell.  Fig¬ 
ure  13  shows  the  four  main  components  of  a  typical  TPV  system. 
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Figure  13.  TPV  System 
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1.  Heat  Source 

The  first  component  is  a  heat  source.  The  temperature  of  the  heat  source  is  a  very 
important  factor  in  the  efficiency  of  the  system.  This  relationship  will  be  discussed  in 
more  depth  in  the  next  section. 

2.  Selective  Emitter 

The  second  component  in  the  system  is  a  selective  emitter.  The  purpose  of  this 
emitter  is  to  help  shape  the  radiation  spectrum  from  the  heat  source.  As  with  solar  cells, 
radiation  absorbed  by  the  cell  that  does  not  have  enough  energy  to  create  an  electron- hole 
pair  just  adds  heat  to  the  cell,  thereby  lowering  the  efficiency.  Rare  earth  oxides  such  as 
ytterbium  oxide  (YblOa)  and  yttrium  aluminum  garnet  (YAG)  have  proven  to  be  very 
effective  materials  because  they  emit  radiation  in  a  line  spectrum  instead  of  a  continuous 
spectrum  like  a  heat  source  but  they  are  currently  in  the  experimental  stages  and  are  not 
always  found  in  TPV  systems  [Ref.  10].  A  selective  emitter  was  not  used  in  these  simu¬ 
lations  because  the  response  of  the  materials  used  was  desired  over  a  wide  range  of  the 
spectrum.  This  data  is  important  when  designing  multijunction  cells  which  will  be  dis¬ 
cussed  later. 


3.  Filter 

The  filter  performs  a  similar  function  to  the  selective  emitter.  It  is  designed  to  re¬ 
flect  incident  radiation  that  does  not  have  enough  energy  to  create  an  electron-hole  pair  in 
the  TPV  cell.  The  radiation  is  reflected  back  to  the  emitter.  This  helps  to  keep  the  tem¬ 
perature  of  emitter  as  high  as  possible  which  improves  the  efficiency  of  the  system  [Ref. 
11]. 

Figure  14  shows  an  experimental  resonant  gold  film  filter  [Ref.  1 1].  Designed  by 
Energetics  Development  Technology  (EDTEK),  this  filter  is  etched  from  a  solid  gold 
film  in  a  pattern  designed  to  resonate  at  a  desired  frequency  [Ref.  11].  The  filter  can  be 
tuned  to  resonate  at  any  desired  frequency  depending  on  the  intended  spectrum.  Tests 
have  shown  that  the  filter  has  very  good  transmission  in  the  desired  band  and  reflects 
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over  98%  of  out-of-band  radiation  [Ref.  11].  Filters  were  not  included  in  the  models 
built  here  for  the  same  reasons  that  selective  emitters  were  not  used. 
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Figure  14.  Resonant  Filter  [After  Ref.  11.] 


Figure  15  shows  the  frequency  response  of  a  filter  built  with  this  technique  along 
with  the  AMO  spectrum  [Ref.  11].  This  particular  filter  was  designed  to  resonate  at  two 
distinct  frequencies.  This  capability  is  useful  for  multijunction  cells  which  consist  of  in¬ 
dividual  cells  that  are  designed  to  respond  to  different  wavelengths  of  light. 


Wavelength  (p.m) 

Figure  15.  Resonant  Filter  Frequency  Response  [From  Ref.  11.] 
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4.  Photo-Electric  Converter 

Conventional  RTG’s  use  a  thermocouple  to  convert  heat  into  electricity  [Ref.  13]. 
A  thermocouple  produces  power  using  what  is  known  as  the  Seebeck  Effect.  Seebeck 
discovered  that  when  two  materials  from  the  thermoelectric  series  were  connected,  a  po¬ 
tential  could  be  measured  across  them  if  the  two  materials  were  at  different  temperatures 
[Ref.  12].  In  addition,  if  the  two  materials  were  connected  with  a  wire  so  as  to  close  the 
circuit,  a  small  current  was  forced  to  flow.  In  an  RTG,  a  SiGe  thermocouple  is  placed  in 
the  vicinity  of  the  heat  source  with  one  side  exposed  to  the  heat  and  the  other  side  to  a 
heat  sink  [Ref.  13].  The  electrical  power  produced  is  proportional  to  the  difference  in 
temperature  between  the  hot  and  cold  junctions.  These  devices  do  a  relatively  poor  job  of 
power  conversion,  however,  only  achieving  efficiencies  of  a  few  percent  in  optimum 
working  conditions  [Ref.  12]. 

This  thesis  will  explore  the  thermophotovoltaic  cell  as  a  candidate  for  replacing 
the  thermocouple  in  RTG’s.  Such  a  cell  operates  in  a  manner  very  similar  to  how  a  solar 
cell  works.  The  difference  between  the  two  is  the  band  of  the  electromagnetic  spectrum 
the  cell  is  designed  to  respond  to.  Where  a  solar  cell  responds  to  the  visible  band  of  the 
spectrum,  thermophotovoltaic  cells  respond  to  the  infrared  region. 

In  order  for  the  TPV  cell  to  respond  to  radiant  energy  from  an  RTG  source,  the 
energy  of  the  radiant  photons  must  be  greater  than  the  energy  bandgap  of  the  material. 

As  previously  mentioned,  photon  energy  is  inversely  proportional  to  the  temperature  of 
the  radiant  body.  From  Table  1  it  can  be  seen  that  the  radiation  from  a  1300K  blackbody 
source  peaks  at  2.23  pm  which  is  a  much  longer  wavelength  than  the  0.45  pm  the  sun’s 
spectrum  peaks  at.  This  means  that  materials  used  for  TPV  cells  must  have  very  low  en¬ 
ergy  bandgaps  in  order  to  respond  to  such  long  wavelengths  of  light.  Table  2  shows  sev¬ 
eral  different  semiconductor  materials  with  their  respective  energy  bandgaps  and  the 
maximum  wavelength  of  light  they  can  respond  to  as  calculated  from  Equation  2.2. 
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Table  2.  Maximum  Frequency  Response  of  Various  Materials 


In  order  for  the  peak  of  the  1300K  blackbody  spectrum  to  occur  in  the  energy 
band  that  a  TPV  cell  can  respond  to,  the  cell  must  be  made  of  a  material  with  an  energy 
bandgap  lower  than  0.55  eV  .  As  can  be  seen  from  Table  2  there  are  materials  that  have 
such  a  low  bandgap;  however  there  are  tradeoffs  to  be  made.  The  open-circuit  voltage  of 
a  cell  is  closely  related  to  the  energy  bandgap  of  the  material.  As  the  energy  bandgap 
goes  down,  so  does  the  open-circuit  voltage.  A  material  with  a  bandgap  of  0.39  eV  may 
only  produce  an  open-circuit  voltage  of  just  a  few  tens  of  millivolts  or  less.  Such  a  low 
voltage  drastically  reduces  the  maximum  power  from  the  cell  in  turn  reducing  the  effi¬ 
ciency. 

Gallium  antimonide  (GaSb)  is  the  most  commonly  published  material  used  for 
TPV  cells.  Although  it  does  not  have  a  low  enough  energy  bandgap  to  respond  to  most 
of  the  radiant  energy  from  an  RTG  heat  source,  the  cells  are  relatively  easy  to  fabricate 
and  are  capable  of  producing  a  significant  amount  of  power  per  square  centimeter.  In¬ 
dium  gallium  arsenide  (InGaAs)  is  another  popular  material  for  TPV  cells.  As  a  ternary 
compound,  many  of  its  electrical  and  optical  properties  including  energy  band  gap 
change  with  varying  molar  concentrations  of  indium  and  gallium.  As  a  0.55  eV  cell,  its 
low  open-circuit  voltage  is  not  very  suitable  for  single  junction  cells,  but  it  is  attractive 
for  use  in  multijunction  cells  [Ref.  14].  For  these  reasons,  GaSb  and  InGaAs  are  the  ma¬ 
terials  that  will  be  studied  in  this  thesis. 
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This  chapter  gave  an  in  depth  discussion  of  thermophotovoltaics  including  the 
cells  themselves  and  the  spectrum  they  are  designed  to  respond.  The  next  chapter  will 
cover  the  systems  that  thermophotovoltaic  cells  are  designed  to  operate  in,  radioisotope 
thermoelectric  generators. 
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IV.  RADIOISIOTOPE  THERMOELECTRIC  GENERATOR 


This  chapter  presents  the  radioisotope  thermoelectric  generator  (RTG).  As  a 
source  of  space  nuclear  power,  the  use  of  RTG’s,  while  proven  safe,  reliable,  and  neces¬ 
sary  for  deep  space  missions,  has  been  greatly  debated.  The  history  of  RTG’s  will  be 
discussed  including  an  account  of  their  successes.  Environmental  issues  regarding  the 
use  of  RTG’s  will  also  be  mentioned  along  with  the  many  measures  in  place  to  maintain 
their  impeccable  safety  record.  Finally  this  chapter  also  includes  about  possible  im¬ 
provements  to  RTG’s  necessary  for  higher  efficiency  operation. 


A.  HISTORY 


As  space  missions  have  begun  to  take  man-made  vehicles  and  satellites  farther 
and  farther  from  the  sun,  the  need  for  a  better  means  of  power  generation  has  become 
more  important.  The  intensity  of  light  from  the  sun  reduces  with  the  square  of  the  dis¬ 
tance  from  it  meaning  that  massive  solar  arrays  would  be  needed  in  order  to  produce 
minimal  power  on  distant  missions.  NASA  has  calculated  that  on  missions  such  as  Cas¬ 
sini,  the  weight  of  solar  arrays  needed  to  provide  the  desired  amount  of  power  would 
more  than  double  the  weight  of  the  entire  satellite  [Ref.  3]. 

An  alternative  to  these  massive  arrays  is  to  bring  a  fuel  supply  that  will  be  capable 
of  providing  adequate  power  over  the  entire  lifespan  of  the  vehicle.  It  is  not  feasible  to 
use  just  batteries  because  the  number  of  batteries  required  for  such  a  long  duration  would 
be  far  too  heavy  [Ref.  3].  Conventional  fuels  are  also  not  viable  options  due  to  the  poor 
mass  to  power  ratios  [Ref.  3]. 

Since  1961,  the  United  States  Department  of  Energy  (DOE)  has  supplied  NASA 
with  RTG’s  powered  by  Pu-238  [Ref.  13].  Table  3  lists  some  missions  that  have  made 
use  of  these  RTG’s  [Ref.  15]. 
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Date 


Spacecraft 


Power 

Source 


Misson  Type 


Status 


i  i 

i  i 

i  i 

29  JUN  1961  j  Transit  4A  |  SNAP-3B7 

i  i 

i  i 

i  i 

:  : 

Navigational 

Operated  for  15  years. 

Shut  down  but  still 

operational. 

. ? . ? . 

i  : 

15  NOV  1961  Transit  4B  SNAP-3B8 

1  1 

i  i 

Navigational 

RTG  operated  for 

9  years 

. ] . \ . 

1  1 

i  : 

i  i 

i  i 

28  SEP  1963  Transit  5-BN-l  SNAP-9A 

1  1 

1  1 

1  ! 

i  i 

i  i 

!  ! 

Navigational 

RTG  operated  as 
planed.  Electrical 
problems  caused 

failure  after  9  months. 

1  1 

1  1 

1  1 

05  DEC  1963  |  Transit  5-BN-2  j  SNAP-9  A 

i  i 

i  i 

i  i 

i  i 

Navigational 

RTG  operated  for  6 
years.  Navigation  failed 
after  1.5  years. 

. \ . \ . 

i  i 

i  i 

21  APR  1964  Transit  5-BN-3  SNAP-9A 

1  1 

1  1 

1  1 

i  i 

Navigational 

Mission  aborted.  RTG 

burned  up  on  reentry 
as  designed. 

i  i 

03  APR  1965  Snapshot  SNAP-10A 

1  1 

1  1 

Experimental 

Successfully  achieved 
orbital  operations. 

. 1 . 1 . 

1  1 

j  j 

18  MAY  1968  I  Nimbus-B-1  !  SNAP-19B2 

:  : 

i  i 

i  i 

i  i 

Meterological 

Mission  aborted. 

RTG  recovered  and 

recycled. 

i  i 

14  APR  1969  !  Nimbus  III  i  SNAP-19B3 

i  i 

i  i 

Meterological 

RTG's  operated  for 
more  than  2.5  years. 

i  i 

14  NOV  1969!  Apollo  12  |  SNAP-27 

!  ! 

Lunar  Surface 

RTG  shut  down  after 

8  years  of  operation. 

i  i 

11  APR  1970  i  Apollo  13  |  SNAP-27 

i  i 

i  j 

Lunar  Surface 

Mission  aborted.  Heat 

Source  returned. 

1  1 

31  JAN  1971  Apollo  14  SNAP-27 

1  1 

:  i 

Lunar  Surface 

RTG  shut  down  after 

6.5  years  of  operation. 

. i . i . 

i  i 

26  JUL  1971  Apollo  15  SNAP-27 

1  1 

i  i 

Lunar  Surface 

RTG  shut  down  after 

6  years  of  operation. 

. ] . ] . 

1  1 

1  1 

02  MAR  1972  Pioneer  10  SNAP-19 

1  1 

1  1 

1  1 

1  1 

Planetary 

RTG  still  operation. 
Spacecraft  now 
beyond  Pluto. 
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i  i  Power  i 

Date  |  Spacecraft  j  j  Mission  Type 

j  Source 

Status 

1  1  1 

16  APR  1972  j  Apollo  16  j  SNAP-27  j  Lunar  Surface 

i  i  i 

i  i  i 

RTG  shut  down  after 

5.5  years  of  operation. 

02  APR  1972  !  Transit  !  Transit-RTG  !  Navigational 

RTG  still  operating. 

"if  f 

07  DEC  1972  Apollo  17  SNAP-27  Lunar  Surface 

i  i  i 

!  !  ! 

RTG  shut  down  after 

5  years  of  operation. 

1  1  1 

05  APR  1973  I  Pioneer  11  !  SNAP- 19  !  Planetary 

ii  i 

i  i  i 

Spacecraft  operated  to 
Jupiter  and  Saturn. 

i  i  i 

20  AUG  1975  i  Viking  1  |  SNAP-19  |  Mars  Surface 

i  i  i 

i  i  i 

RTG  shut  down  after 

6  years  of  operation. 

i  i  i 

09  SEP  1975  |  Viking  2  j  SNAP- 19  j  Mars  Surface 

i  1  1 

1  1  1 

RTG  operated  4  years 
until  relay  link  was  lost. 

14  MAR  1976  LES  8  MHW-RTG  Communication 

!  !  ! 

RTG  still  operating. 

14  MAR  1976  LES  9  MHW-RTG  Communication 

RTG  still  operating. 

. 1 . 1 . 1 . 

1  1  1 

j  j  j 

20  AUG  1977  Voyager  2  MHW-RTG  Planetary 

i  1  1 

1  1  1 

1  1  1 

RTG  still  operating. 
Operated  to  Jupiter, 
Saturn,  Uranus,  Neptune. 

i  i  i 

i  i  i 

i  i  i 

05  SEP  1977  Voyager  1  MHW-RTG  Planetary 

1  1  1 

i  i  i 

i  i  i 

!  !  ! 

RTG  still  operating. 
Operated  to  Jupiter, 
Saturn,  and  beyond. 

1  1  1 

1  1  1 

1  1  1 

18  OCT  1989  i  Galileo  |  GPHS-RTG  |  Planetary 

i  i  i 

i  i  i 

i  i  i 

i  i  i 

RTG  still  operating. 
Spacecraft  orbiting 
around  Jupiter. 

. ? . \ . \ . 

i  i  i 

i  i  i 

06  OCT  1990  |  Ulysses  j  GPHS-RTG  Planetary/Solar 

i  i  i 

i  i  i 

i  i  i 

RTG  still  operating. 
Spacecraft  in  polar 

orbit  around  the  sun. 

i  i  i 

|  |  | 

15  OCT  1997  j  Cassini  j  3-RTGs  j  Planetary 

1  1  1 

1  1  1 

i  i  i 

i  i  i 

Table  3.  RTG  Missions  [After  Ref.  15. 

RTGs  still  operatin. 
Spacecraft  en  route 

to  Saturn. 

] 

This  table  shows  that  RTG’s  are  capable  of  providing  power  over  long  periods  of 
time.  The  first  RTG  sent  into  space  in  June  of  1961  on  the  Transit  4A  spacecraft  is  still 
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providing  power  [Ref.  15].  Figure  16  shows  a  rendering  of  a  current  RTG  module  [Ref. 
16]. 
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Figure  16.  RTG  Module  [From  Ref.  16.] 


B.  FUEL  SUPPLY 


An  RTG  fuel  supply  must  be  able  to  balance  a  high  power- to- mass  ratio  with  a 
long  lifetime  and  reasonable  safety.  As  mentioned  previously,  current  RTGs  have  used 
Pu-238  as  a  heat  source.  Pu-238  has  a  half-hfe  of  87  years,  meaning  that  after  20  years  of 
operation  it  will  still  produce  80%  of  the  power  it  produced  when  new  [Ref.  17].  Figure 
17  shows  the  power  output  versus  time  for  a  Pu-238  heat  source  [Ref.  17].  This  isotope, 
when  pure,  has  a  relatively  large  power  to  mass  ratio  of  540  W/kg  [Ref.  17]. 


32 


i  Years 

Figure  17.  Pu-238  Power  Output  [From  Ref.  17.] 

Other  fuel  supplies  have  been  considered  for  use  with  RTG’s  as  well.  Strontium 
90  (Sr- 90)  for  example  was  considered  as  a  possible  fuel  source.  With  a  half  life  of  28. 1 
years,  Sr- 90  decays  much  faster  than  Pu-238  [Ref.  18].  After  20  years,  it  will  only  pro¬ 
duce  approximately  60%  of  the  power  it  originally  did  [Ref.  18].  In  addition,  the  power 
to  mass  ratio  when  new  is  only  460  W/kg  which  is  more  than  15%  lower  that  that  of  Pu- 
238  [Ref.  18], 

Although  the  use  of  Pu-238  has  startled  many  people,  it  is  actually  a  rather  safe 
fuel.  In  order  to  demonstrate  to  the  public  how  safe  the  use  of  Pu-238  is,  President  Ei¬ 
senhower  had  the  picture  in  Figure  18  published  in  the  newspaper  taken  in  the  Oval  Of¬ 
fice  at  the  White  House  with  senior  members  United  States  Atomic  Energy  Commission 
(AEC)  and  an  early  RTG  module  [Ref.  17]. 
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C.  SAFETY 


1.  Accident  History 

In  total  the  United  States  Department  of  Energy  has  provided  44  RTG’s  and  more 
than  240  general-purpose  heat  source  (GPHS)  units  since  1961  with  only  three  accidents 
[Ref.  17].  The  first  of  these  accidents  occurred  on  the  Transit  5-BN-3  launched  in  April 
of  1964.  This  mission  had  to  be  aborted  when  the  launch  vehicle  failed.  The  vehicle  and 
RTG  burned  up  in  the  atmosphere  upon  reentry  where  the  Pu-238  fuel  source  was  vapor¬ 
ized  and  dispersed  into  the  upper  atmosphere  as  it  was  designed  to  do  [Ref.  17].  Shortly 
after  this  accident,  the  RTG  module  was  redesigned  to  survive  atmospheric  reentry  with 
the  fuel  source  intact  [Ref  .17]. 

The  second  accident  occurred  on  the  Nimbus -B-l  launched  in  May  of  1968  [Ref. 
17].  This  mission  had  to  be  aborted  shortly  after  launch  by  a  range  safety  destruct  [Ref. 
17].  The  newly  designed  RTG  module  held  up  through  reentry  and  was  recovered  with¬ 
out  any  radiation  leakage.  The  fuel  source  was  reused  on  later  missions  [Ref.  17]. 

The  third  and  final  accident  occurred  on  Apollo  13  launched  in  April  of  1970. 
This  mission  carried  an  RTG  intended  to  provide  power  to  a  lunar  seismic  station  [Ref. 
17].  The  RTG  was  located  on  the  lunar  module  which  broke  up  on  reentry  into  earth’s 
atmosphere  [Ref.  17].  The  RTG  remained  intact  and  is  now  located  at  the  bottom  of  the 
Tonga  Trench  in  the  Pacific  Ocean  [Ref.  17]. 
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2.  Safety  Measures 

Many  features  have  been  added  to  RTG’s  to  ensure  their  safe  operation  under  all 
conditions.  Figure  19  shows  a  diagram  of  a  modem  RTG  unit  [Ref.  19]. 
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Figure  19.  RTG  Safety  Measures  [From  Ref.  19.] 


It  can  be  seen  from  the  figure  that  the  fuel  pellet  makes  up  a  very  small  percent¬ 
age  of  the  unit  while  added  safety  features  occupy  the  majority  of  the  area.  The  fuel  itself 
has  even  been  modified  for  safety  purposes.  The  original  Pu-238  metal  has  been  replaced 
with  plutonium  dioxide  [Ref.  17].  This  heat  resistant  ceramic  form  prevents  the  fuel  from 
vaporizing  under  extreme  heat  conditions  such  atmospheric  reentry  [Ref.  17].  The  mate¬ 
rial  is  also  insoluble  and  has  a  very  low  chemical  reactivity  making  it  safe  even  if  it  ends 
up  in  the  ocean  like  the  module  from  Apollo  13  [Ref.  17].  This  ceramic  form  also  influ¬ 
ences  how  the  material  breaks  up.  Instead  of  becoming  a  fine  powder  that  can  easily 
spread,  it  fractures  into  large  non- respirable  chunks  [Ref.  17].  As  an  added  safety  meas¬ 
ure,  the  fuel  is  separated  into  small  pellets  that  are  each  individually  encapsulated  in  a 
heat  shield  and  impact  shell.  The  module  itself  is  enclosed  in  an  iridium  capsule  and  sur¬ 
rounded  by  graphite  blocks  [Ref.  17] .  These  are  high  strength  and  corrosion  and  heat  re¬ 
sistant  materials.  The  RTG  unit  is  designed  to  withstand  conditions  it  could  never  realis¬ 
tically  encounter. 
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The  Pu-238  in  RTG’s  is  an  alpha  emitter  [Ref.  17].  The  majority  of  alpha  parti¬ 
cles  can  be  absorbed  by  a  piece  of  paper  or  the  epidermal  layer  of  human  skin  [Ref.  17]. 
RTG’s  do  not  utilize  either  fission  of  fusion  in  producing  heat  [Ref.  17].  Because  of  this, 
an  RTG  could  never  explode  like  a  nuclear  bomb.  In  addition,  no  accident  involving  an 
RTG  could  ever  create  the  acute  radiation  sickness  associated  with  nuclear  explosions  or 
meltdowns  [Ref.  17]. 

3.  Safety  Testing 

As  previously  mentioned,  RTG  modules  are  designed  to  endure  conditions  they 
could  never  realistically  face.  In  order  to  ensure  the  safety  of  the  modules,  they  are  tested 
under  extremely  harsh  conditions  [Ref.  17]. 


a.  Fire 

Extreme  heat  upon  atmospheric  reentry  and  other  unforeseen  accidents 
can  cause  fires  onboard  a  spacecraft.  RTG’s  have  been  directly  exposed  to  intense  fires 
with  minimal  damage  [Ref.  17].  Most  importantly,  no  fuel  was  released  during  exposure. 

b.  Blast 

A  shuttle  explosion  can  cause  tremendous  shockwaves.  The  RTG  module 
was  tested  under  blast  waves  stronger  than  those  expected  from  such  an  explosion  and 
once  again  remained  intact  without  releasing  any  nuclear  fuel  or  harmful  radiation  from 
the  unit  [Ref.  17]. 


c.  Reentry 

In  the  event  that  a  satellite  degrades  in  orbit  enough  to  reenter  the  earths 
atmosphere,  it  will  experience  intense  amounts  of  heat.  An  RTG  is  designed  to  survive 
reentry  while  preventing  the  Pu-238  fuel  source  from  vaporizing  in  the  atmosphere.  This 
environment  was  tested  using  an  arc-jet  furnace,  and  the  RTG  was  able  withstand  the 
heat  safely  [Ref.  17]. 
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d.  Earth  Impact 

Since  RTG’s  were  redesigned  to  survive  atmospheric  reentry  intact,  they 
also  had  to  be  strong  enough  to  survive  impact  with  earth’s  surface.  RTG’s  were  tested 
at  the  maximum  estimated  velocity  of  the  module  approaching  earth’s  surface  of  120 
miles  per  hour.  These  tests  showed  that  in  an  impact  with  water,  soil,  and  sand  the  mod¬ 
ule  would  remain  safely  intact.  An  impact  with  rock  or  concrete,  however,  sometimes 
caused  slight  leakage  resulting  in  small  amounts  of  low-level  radiation  in  localized  areas 
[Ref.  17]. 


e.  Immersion  in  Water 

Because  the  vast  majority  of  earth’s  surface  is  covered  in  water,  RTG’s 
must  be  able  to  survive  in  an  aquatic  atmosphere.  Long  term  exposure  to  sea  water 
proved  the  protective  iridium  capsule  to  be  corrosion  resistant  and  the  fuel  to  be  highly 
insoluble  [Ref.  17]. 

/.  Shrapnel 

In  the  event  of  an  explosion  or  impact  with  space  debris,  an  RTG  could  be 
subjected  to  shrapnel.  Aluminum  and  titanium  bullets  were  used  to  simulate  such  shrap¬ 
nel  from  a  launch  explosion.  Fragments  at  speeds  much  greater  than  those  predicted  to 
be  present  in  such  an  event  were  unable  to  penetrate  the  module  [Ref.  17]. 


g.  Large  Fragments 

In  the  event  of  a  solid  rocket  booster  failure,  the  RTG  module  could  be 
subject  to  impact  by  large  fragments  of  material.  To  simulate  this  scenario,  large  steel 
plates  were  fired  at  an  RTG.  The  studies  showed  that  damage  only  occurred  in  the 
unlikely  event  that  the  comer  or  sharp  edge  of  a  large  fragment  struck  the  module  [Ref. 
17]. 
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D.  FUTURE  IMPROVEMENTS 


Current  RTG’s  use  silicon/germanium  thermocouples  discussed  in  the  previous 
chapter  to  convert  thermal  power  from  the  Pu-238  heat  source  into  electrical  power  [Ref. 
13].  These  thermocouples  operate  at  very  low  efficiency  levels  of  only  a  few  percent. 

The  key  to  improving  RTG  power  production  is  to  find  a  more  efficient  way  to  convert 
thermal  power  into  electrical  power. 

1.  Dynamic  Isotope  Power  Systems 

Dynamic  isotope  power  systems  (DIPS)  convert  thermal  energy  into  electrical  en¬ 
ergy  using  moving  parts.  Created  in  1816  by  Robert  Stirling,  the  Stirling  engine  is  the 
most  important  dynamic  isotope  power  system  to  space  power  systems  [Ref.  20].  Unlike 
combustion  engines,  Stirling  engines  have  a  closed  gas  chamber  [Ref.  20].  There  are  no 
explosions  that  take  place  in  a  Stirling  engine  and  there  is  no  gas  exhausted  from  the  en¬ 
gine.  The  power  to  run  this  engine  comes  from  an  external  heat  source  [Ref.  20].  Figure 
20  shows  a  simplified  example  of  a  Stirling  engine  [Ref.  20]. 


There  are  four  steps  in  the  Stirling  cycle  [Ref.  20] : 

1.  Heat  is  added  to  the  hot  cylinder  (left),  causing  pressure  to  build.  This  pushes 
the  left  piston  down. 
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2.  The  left  piston  moves  up  and  the  right  piston  moves  down.  This  pushes  the 
heated  gas  from  the  left  cylinder  into  the  cold  cylinder  on  the  right  where  the  gas  is 
quickly  cooled  down.  When  the  gas  cools,  the  pressure  goes  down  making  it  easier  to 
compress  in  step  3. 

3.  The  piston  in  the  cold  cylinder  on  the  right  compresses  the  gas.  Heat  gener¬ 
ated  by  this  compression  is  almost  instantly  removed  by  the  cooled  cylinder. 

4.  The  right  piston  moves  up  while  the  left  piston  moves  down  forcing  the  gas 
into  the  heated  cylinder  on  the  right.  This  cylinder  quickly  heats  the  gas  raising  the  pres¬ 
sure  in  the  cylinder. 

The  cycle  then  repeats. 


In  this  four- step  cycle,  power  is  only  created  in  the  first  step.  There  are  two  ways 
to  increase  the  power  output  of  a  Stirling  engine.  The  first  method  is  to  increase  the 
power  output  of  stage  one  [Ref.  20].  This  is  accomplished  by  raising  the  pressure  of  the 
gas  in  the  hot  cylinder  which  is  done  by  increasing  the  temperature  of  the  external  heat 
source.  The  second  method  for  increasing  the  power  output  is  lowering  the  gas  pressure 
in  the  cold  cylinder  [Ref.  20].  This  is  accomplished  by  lowering  the  temperature  of  the 
cylinder.  The  Stirling  engine,  like  a  thermocouple,  produces  a  power  output  related  to 
the  difference  in  temperature  of  the  hot  and  cold  cylinders. 

The  Stirling  engine  is  not  suitable  for  widespread  use  for  two  main  reasons.  Be¬ 
cause  the  heat  source  is  external  to  the  engine,  it  takes  some  time  for  the  engine  to  warm 
up  before  it  can  produce  useful  power  [Ref.  20].  Secondly,  the  engine  is  slow  to  respond 
to  changes  in  the  power  supply  [Ref.  20]. 

Stirling  engines  do  however,  show  promise  as  a  means  of  power  generation  for 
space  based  systems.  NASA  Glenn  has  been  developing  Stirling  engine  technology  since 
the  1970’s  and  has  been  developing  radioisotope  power  systems  for  deep  space  missions 
since  1990  [Ref.  21].  Figure  21  shows  a  Stirling  generator  developed  by  NASA  Glenn 
[Ref.  21]. 
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Figure  21.  110  Watt  Electric  Stirling  Generator  [From  Ref.  21.] 

This  figure  shows  two  synchronized  55  watt  electric  generators.  This  configura¬ 
tion  reduces  vibrations  in  the  system.  The  addition  of  internal  precision  springs  and 
dampers  causes  this  engine  to  operate  at  a  resonant  frequency  that  increases  the  effi¬ 
ciency  of  the  system.  This  model  has  been  tested  to  work  at  an  efficiency  level  of  over 
20%  [Ref.  21]. 

2.  Alkaline  Metal  Thermal -to- Electric  Conversion 

Figure  22  shows  the  schematic  diagram  of  a  simple  alkaline  metal  thermal-to- 
electric  conversion  (AMTEC)  system  [Ref.  22]. 


A 


B 


C 


A- Heat  Source  B-  Liquid  Sodium  C-  Heat  Sink 
D-  Sodium  Vapor  E-  Beta-Alumina  Solid  Electrolyte 


Figure  22.  AMTEC  Schematic  [From  Ref.  22.] 
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The  key  to  this  converter  is  the  beta-alumina  solid  electrolyte  (BASE).  This  con¬ 
ducts  positive  sodium  ions  much  better  than  it  conducts  sodium  atoms  or  electrons  [Ref. 
22].  When  heat  is  applied  to  one  side,  a  sodium  pressure  difference  occurs  across  the 
base  and  the  positive  ions  move  through  the  base  and  collect  on  the  low  pressure  side 
leaving  electrons  to  collect  on  the  high  pressure  side  [Ref.  22].  This  creates  a  potential 
across  the  base  that  can  force  a  current  to  flow  through  a  load  [Ref.  22].  Similar  to  the 
previously  discussed  thermoelectric  converters,  the  power  produced  by  an  AMTEC  sys¬ 
tem  is  dependent  on  a  large  temperature  difference  between  the  hot  and  cold  ends,  but 
currently  the  efficiency  of  such  systems  is  very  low  [Ref.  22]. 

3.  Thermophotovoltaic  Cells 

Thermophotovoltaic  (TPV)  cells  are  currently  the  most  promising  thermoelectric 
converters  for  replacing  the  thermocouple  in  RTG’s.  Although  they  are  currently  not  ca¬ 
pable  of  the  operating  efficiencies  of  Stirling  engines,  TPV  cells  have  the  advantage  of  no 
moving  parts  [Ref.  3].  These  cells  are  discussed  in  more  depth  in  the  previous  chapter. 

The  radioisotope  thermoelectric  generator  was  introduced  in  this  chapter.  The 
current  system  as  well  as  possible  future  improvements  were  covered.  The  next  chapter 
will  begin  the  discussion  on  the  simulation  software  used  in  this  thesis. 
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V.  SIMULATION  SOFTWARE 


There  is  currently  a  great  deal  of  research  being  done  on  building  accurate  models 
of  solar  and  thermophoto  voltaic  cells.  While  some  models  can  provide  relatively  accu¬ 
rate  results  in  limited  cases,  there  is  no  widely  available  software  package  that  is  capable 
of  modeling  all  of  the  necessary  parameters  needed  to  provide  models  of  cells  of  various 
materials  under  many  different  spectra.  This  thesis  explores  the  use  of  Silvaco  Virtual 
Wafer  Fabrication  software  to  accomplish  this  task. 


A.  SILVACO 


Silvaco  International  is  a  software  company  that  provides  programs  for  modeling 
all  areas  of  electronics  including  digital  and  analog  circuits  [Ref.  23].  The  company  has 
software  ranging  from  simple  Spice  modeling  to  advanced  integrated  circuit  layout  and 
extraction  tools.  The  Virtual  Wafer  Fabrication  package  is  part  of  Silvaco’s  TCAD  suite. 
It  was  designed  to  accurately  emulate  physical  wafer  manufacturing  [Ref.  23].  The  pack¬ 
age  allows  users  to  easily  optimize  all  aspects  of  IC  design  before  actually  fabricating 
them.  ATLAS  is  the  main  program  in  this  Virtual  Wafer  Fabrication  package  that  will  be 
used  in  this  research.  This  program  can  handle  one-,  two-,  and  three-dimensional  designs 
and  contains  the  following  comprehensive  models  [Ref.  23] : 

•  DC,  AC  small-  signal,  and  full  time-dependency 

•  Drift- diffusion  transport  models 

•  Energy  balance  and  hydrodynamic  transport  models 

•  Lattice  heating  and  heatsinks 

•  Graded  and  abrupt  heterojunctions 

•  Optoelectronic  interactions  with  general  ray  tracing 

•  Amorphous  and  polycrystalline  materials 


General  circuit  environments 
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Stimulated  emission  and  radiation 


•  Fermi-Dirac  and  Boltzmann  statistics 

•  Advanced  mobility  models 

•  Heavy  doping  effects 

•  Full  acceptor  and  donor  trap  dynamics 

•  Ohmic,  Schottky,  and  insulating  contacts 

•  SRH,  radiative.  Auger,  and  surface  recombination 

•  Impact  ionization  (local  and  non-  local) 

•  Floating  gates 

•  Band-to-band  and  Fowler-Nordheim  tunneling 

•  Hot  carrier  injection 

•  Thermionic  emission  currents 

ATLAS  is  a  powerful  program,  but  it  relies  on  the  aid  of  several  other  Silvaco 
programs  to  build  an  accurate  model  of  a  circuit.  The  program  Deckbuild  is  the  run  time 
environment  for  this  software  package.  A  TPV  cell  can  be  defined  graphically  from  the 
command  line  in  this  program  in  a  method  that  will  be  discussed  later.  ATHENA,  the 
process  simulator,  actually  builds  the  devices  defined  by  the  commands  entered  in  the 
Deckbuild  interface  [Ref.  23].  DevEdit  is  another  program  that  allows  a  user  to  define 
cell  graphically  by  defining  mesh  lines  and  nodes  on  a  plane  [Ref.  23].  This  can  also  be 
done  directly  through  Deckbuild  as  will  be  discussed  later.  In  addition  to  run  time  out¬ 
put,  ATLAS  also  writes  *.log  files  which  contain  the  output  data  solved  for  during  simu¬ 
lations.  The  output  gathered  from  these  simulations  was  voltage  and  current  levels  and 
photogeneration  rates.  The  included  program  TonyPlot  can  input  these  *.log  files  and 
present  the  information  in  a  graphical  form.  Figure  23  shows  a  block  diagram  of  how 
these  programs  interact  to  solve  a  particular  TPV  cell  [Ref.  23] . 
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Figure  23.  Simulator  Block  Diagram  [From  Ref.  23.] 


B.  DEFINING  A  CELL 


In  this  thesis,  solar  and  TPV  cells  were  defined  graphically  in  ATLAS  through  the 
Deckbuild  interface.  In  order  to  fully  define  a  cell,  it  was  necessary  to  enter  not  just  di¬ 
mensions,  but  also  material  properties.  The  following  format  was  used  in  building  each 
model. 

•  Defining  constants  (cell  dimensions,  doping  levels,  shading,  etc.) 

•  Defining  mesh  boundaries  using  the  previously  defined  constants 

•  Building  the  X-  mesh 

•  Building  the  Y- mesh 

•  Defining  regions  (anode,  cathode,  contacts,  etc.) 

•  Placing  electrodes 

•  Defining  doping  levels  in  each  region 

•  Defining  material  properties  (electrical,  optical,  band,  etc.) 

•  Adding  incident  spectrum 

•  Solving 
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1.  Defining  Constants 

ATLAS  allows  a  user  the  option  of  building  a  model  directly  by  inputting  dimen¬ 
sions  in  numbers  or  indirectly  by  entering  dimensions  as  variables  names  that  have  been 
previously  defined.  In  this  thesis,  it  was  necessary  to  be  able  to  change  various  dimen¬ 
sions  quickly  and  easily  so  the  second  method  of  defining  a  cell  was  utilized. 

The  first  portion  of  the  code  was  dedicated  to  defining  constants  that  would  later 
be  used  in  calculations.  For  example 

set  cathodethickness=10 

set  cathodedoping=1.000000e+020 

set  anodethickness=.37 

set  anodedoping=9.8el9 

set  bsfthickness=  15. 000000 

set  bsfdoping=l  .000000e+020 

defines  constants  that  represent  the  thickness  and  doping  level  for  the  cathode,  anode,  and 
BSF  layers  of  a  cell.  These  variables  can  now  be  used  later  in  the  program.  The  benefit 
to  defining  dimensions  in  this  manner  is  that  by  simply  changing  the  variable  definition, 
the  whole  model  is  altered  to  reflect  this  change  making  iterations  quick  and  easy  to  exe¬ 
cute. 


2.  Defining  Mesh  Boundaries 

ATLAS  solves  equations  at  many  different  nodes  in  a  defined  area.  These  nodes 
are  defined  by  meshes.  The  constants  defined  in  the  first  part  of  the  ATLAS  code  are 
used  in  defining  mesh  boundaries.  The  following  commands  are  an  example  of  how 
mesh  boundaries  are  defined  in  a  cell. 


set  bsfL=0 

set  bsfH=$bsfL-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 
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set  catL=$bsfH 


set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/$divs 
set  anoL=$catH 

set  anoH=$anoL-$anodethickness 
set  anodiv  s=$anodethicknes  s/$div  s 

The  dollar  sign  is  an  indication  to  ATLAS  that  the  following  string  is  a  user  de¬ 
fined  constant.  In  ATLAS,  moving  up  on  the  Y-axis  is  in  the  negative  direction  so  these 
commands  define  a  cell  by  starting  at  Y  =  0  and  building  upwards  including  divisions  in 
each  region. 


3.  Building  the  X-  and  Y-  Meshes 

ATLAS  requires  that  each  region  boundary  fall  on  a  mesh.  In  order  to  ensure  that 
this  is  the  case,  mesh  lines  are  manually  drawn  at  each  of  these  boundaries  using  the  fol¬ 
lowing  commands. 

x.mesh  loc=0  spac=$cellwdiv 

x.mesh  loc=$contLL  spac=$capwdiv 

x.mesh  loc=$contRR  spac=$capwdiv 

x. mesh  loc=$cellwidth  spac=$cellwdiv 

y. mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 
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Each  of  these  command  lines  defines  both  a  location  and  a  spacing  interval. 
ATLAS  will  draw  mesh  lines  at  each  defined  location  as  well  as  lines  at  the  defined  spac¬ 
ing  interval.  When  two  defined  mesh  lines  with  different  spacing  intervals  are  next  to 
each  other,  ATLAS  will  automatically  grade  the  mesh  spacing  between  the  two  lines. 
Ligure  24  shows  an  example  of  a  typical  mesh. 


Ligure  24.  Typical  Mesh 


When  defining  a  mesh,  it  is  important  to  take  careful  note  of  the  spacing.  In  order 
to  get  an  accurate  solution,  the  mesh  must  be  fine  enough  to  get  an  accurate  representa¬ 
tion  of  the  cell.  If  the  mesh  is  too  fine,  however,  it  will  take  an  unnecessarily  long  time  to 
solve  all  of  the  nodes. 

4.  Defining  Regions 

Each  layer  of  a  cell  must  be  defined  by  a  material  and  region  boundaries.  As  pre¬ 
viously  mentioned,  these  boundaries  must  fall  on  already  defined  mesh  lines.  The  fol¬ 
lowing  commands  arc  used  to  define  regions  in  the  GaSb  cell  shown  in  Ligure  25.  The 
bracketed  numbers  represent  the  value  of  each  variable  used  in  micrometers  and  are  not 
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included  in  the  actual  code.  This  cell  has  a  0.1  pm  thick  cap,  0.37  pm  anode,  and  a 
10  pm  thick  cathode.  The  BSF  layer  defined  in  these  commands  cannot  be  seen  in  Fig¬ 
ure  25. 


#Anode 

region  num=5  material=GaSb  x.min=0  x.max=$cellwidth[500] 
y.min=$anoH[-350.37]  y.max=$anoL[-350] 

#Cathode 

region  num=6  material=GaSb  x.min=0  x.max=$cellwidth[500] 
y  .min=$catH  [-350]  y  .max=$c  atL  [-340] 

#BSF 

region  num=7  material=GaSb  x.min=0  x.max=$cellwidth[500] 
y.min=$bsfH[-340]  y.max=$bsfL[0] 

Figure  25  shows  a  plot  of  cell  regions  generated  by  TonyPlot. 
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Figure  25.  Cell  Regions 
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5.  Placing  Electrodes 

Once  the  cell  regions  are  fully  defined,  electrodes  must  be  defined  to  electrically 
connect  to  the  cell.  ATLAS  allows  a  user  to  define  an  ideal  contact  or  specify  a  partic  u- 
lar  material  with  non- ideal  properties.  It  is  even  possible  to  add  external  discrete  compo¬ 
nents  to  the  cell  such  as  resistors,  inductors,  and  capacitors.  Electrodes  are  defined  using 
the  following  commands. 

electrode  name=anode  material=Gold  x.min=$contLL  x.max=$contRR 
y.min=$contH  y.max=$contL 

electrode  name=cathode  material=Gold  x.min=0  x.max=$cellwidth 
y.min=$bsfL  y.max=$bsfL 

6.  Defining  Doping  Levels 

The  doping  level  in  each  region  can  be  defined  independently  using  the  following 
commands. 

#  Anode 

doping  uniform  region=5  p.type  conc=$anodedoping 

#  Cathode 

doping  uniform  region=6  n.type  conc=$cathodedoping 

This  defines  a  constant  doping  level  in  a  region.  The  boundaries  between  two  re¬ 
gions  will  have  abrupt  junctions.  For  more  realistic  and  more  complex  doping  schemes, 
Silvaco  provides  a  C- Interpreter,  which  allows  a  user  to  define  any  conceivable  doping. 

7.  Defining  Material  Properties 

Using  the  following  commands,  the  user  can  define  as  many  material  properties 
as  are  available. 

material  material=GaSb  COPT=5e-10 

material  material=GaSb  EG300=.726  PERMITTIVITY=14.4 
AFFINITY=4.06 
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material  material=GaSb  NC300=2.1el7  NV300=1.8el9 


material  material=GaSb  index.file=GaSb.opt 

These  commands  define  optical  recombination,  energy  band  gap,  permittivity, 
electron  affinity,  available  conduction  and  valence  band  states,  and  optical  parameters  for 
the  material  gallium  antimonide  (GaSb).  The  index  file  ‘GaSb.opt’  contains  refractive 
index  and  extinction  parameters  for  a  span  of  wavelengths  in  the  area  of  interest.  This  is 
one  of  the  most  important  material  properties  to  accurately  define  and  subsequently  one 
of  the  most  difficult  to  find.  Figure  26  shows  a  plot  of  the  data  from  ‘GaSb.opt’. 


8.  Incident  Spectrum 

ATLAS  allows  a  user  to  define  a  spectrum  with  a  *.spec  file  as  follows. 

beam  num=l  x.origin=0  y.origin=$lightY  angle=90  back.refl 

power.file=1300BB.spec  wavel.start=0.207  wavel.end=5 
wavel.num=1000 
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The  *.spec  file  gives  the  spectral  radiance  at  various  wavelengths  in  the  region  of 
interest.  This  command  tells  ATLAS  to  read  in  the  data  from  the  file  ‘1300BB.spec’  and 
to  solve  a  cell  using  1000  evenly  spaced  rays  in  the  interval  from  0.207  pm  -5  pm  .  The 

command  also  defines  the  origin  of  the  beams  and  the  incident  angle. 

9.  Solving 

The  cell  and  its  environment  are  now  completely  defined  in  ATLAS  and  the  pro¬ 
gram  can  now  use  this  model  to  solve.  The  equations  that  are  to  be  solved  depend  on 
what  type  of  output  is  desired.  This  thesis  concentrates  on  current -voltage  relationships 
and  photogeneration  rates.  By  providing  ATLAS  with  a  number  of  current  levels,  it  can 
solve  for  the  associated  voltage  levels  in  order  to  produce  an  I-V  curve  as  follows. 

solve  ianode=0  bl=$I 

solve  ianode=-$il  b  1  =$1 

solve  ianode=-$i2  b  1  =$1 

solve  ianode=-$i3  bl=$I 

solve  ianode=-$i4  bl=$I 


If  photogeneration  rates  are  desired,  the  user  can  tell  ATLAS  to  solve  for  the  pho- 
togeneration  with  discrete  wavelengths  of  light  incident  on  the  cell  as  follows. 

solve  b  1=0.1  lambda=0.3 

solve  bl=0.1  lambda=0.4 

solve  b  1=0.1  lambda=0.5 

solve  b  1=0.1  lambda=0.6 


The  output  of  the  data  generated  from  solving  these  equation  is  stored  in  a  *.log 
file  which  can  be  viewed  graphically  using  Tony  Plot. 
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10.  Source  Code 

The  sample  code  given  in  this  section  was  taken  from  a  model  for  a  GaSb  cell. 
The  full  source  code  for  both  this  cell  and  an  InGaAs  cell  can  be  found  in  Appendix  A. 


C.  MATLAB 


The  Silvaco  software  package  is  capable  of  performing  all  of  the  actions  neces¬ 
sary  to  model  a  cell  and  graph  the  output,  but  it  was  often  useful  to  use  Matlab  programs 
to  supplement  Silvaco  for  some  operations. 

1.  Iterations 

While  ATLAS  is  capable  of  solving  many  very  complex  equations,  it  does  lack 
some  very  basic  programming  capabilities  such  as  conditional  statements  and  loops  that 
are  necessary  for  performing  iterations.  Several  programs  (found  in  Appendix  B)  were 
written  in  Matlab  to  execute  these  functions.  The  program  ‘atlasrun.m’  was  used  to  run 
the  iterations.  In  each  iteration,  ‘atlasrun.m’  would  run  an  input  deck  with  ATLAS,  read 
and  analyze  the  output  data,  iterate  a  parameter  in  the  input  deck,  and  rerun  ATLAS.  The 
program  ‘ivmaxp.m’  solved  for  open-circuit  voltage,  short  circuit  current,  maximum 
power,  fill  factor,  and  efficiency  for  each  iteration  and  recorded  this  data  in  an  output  file, 
‘filerw.m’  was  used  to  iterate  a  specified  parameter  in  the  input  deck  after  each  run. 

2.  Plotting 

The  TonyPlot  program  provided  by  Silvaco  is  capable  of  producing  plots  of  all 
data  that  is  generated  by  ATLAS,  but  it  often  more  convenient  to  use  Matlab  to  plot  the 
data.  In  addition  to  running  the  iteration  scheme,  ‘atlasrun.m’  would  also  update  an  effi¬ 
ciency  plot  after  each  run  of  ATLAS.  Since  most  optimization  schemes  took  several 
hours  to  run,  this  feature  provided  up  to  date  visual  feedback  on  the  progress  of  the  rou¬ 
tine.  The  program  ‘filerw.m’  was  also  used  to  plot  I-V  curves  with  a  marked  maximum 
power  point  on  each  curve. 

All  of  these  programs  can  be  found  in  full  in  Appendix  B. 
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This  chapter  covered  the  simulation  software  used  in  this  thesis.  Silvaco  and  Mat- 
lab  were  both  utilized  extensively  for  building  and  analyzing  TPV  cell  models.  The  next 
chapter  will  cover  the  actual  models  built  along  with  their  resulting  outputs. 
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VL  MODELS  AND  RESULTS 


This  chapter  discusses  the  models  that  were  built  in  ATLAS.  In  each  case,  a  sam¬ 
ple  model  was  built  and  the  output  of  this  model  was  compared  to  accepted  values.  Once 
the  model  was  verified,  it  was  optimized  for  both  AMO  and  the  spectrum  from  a  1300K 
blackbody  using  an  iteration  routine  written  in  Matlab. 


A.  MODEL  CONFIRMATION 


Before  a  TPV  cell  could  be  optimized  for  a  blackbody  spectrum,  it  was  first  nec¬ 
essary  to  confirm  that  the  model  being  used  was  accurate.  Although  there  was  no  data 
available  on  cells  under  a  1300K  blackbody  spectrum,  there  was  available  data  on  many 
different  cells  under  the  AMO  spectrum.  In  addition,  using  Equation  2.2  the  frequency 
response  of  a  cell  could  be  estimated.  The  assumption  is  that  if  the  model  can  accurately 
predict  the  output  of  a  cell  under  AMO  and  has  the  proper  frequency  response,  it  should 
also  accurately  predict  the  output  under  a  different  spectrum. 

1.  Gallium  Arsenide 

Gallium  Arsenide  (GaAs),  while  not  suitable  for  TPV  applications,  is  one  of  the 
most  widely  used  materials  for  solar  cells  and  therefore  the  maximum  current-voltage 

relationship  for  this  cell  is  very  well  documented.  A  1  cm2  GaAs  cell  should  provide  ap¬ 
proximately  Isc  =31  mA  and  Voc  =  1  V  [Ref.  3].  The  cell  in  Figure  27  was  built  and 

tested  in  ATLAS.  ATLAS  was  found  to  do  a  poor  job  of  modeling  antireflective  coat¬ 
ings,  so  a  5%  loss  was  assumed  in  all  models.  With  present  AR  coating  technology  this 
is  a  very  reasonable  assumption. 


55 


Figure  27.  GaAs  Solar  Cell 


Figure  28  shows  the  IV  curve  from  this  cell. 


I-V  Curve 
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From  Equation  2.2,  the  spectral  response  of  a  GaAs  cell  should  approach  zero  at  a 
wavelength  of  approximately  0.89  pm  .  Figure  29  shows  the  calculated  spectral  response 
with  a  peak  at  0.87  pm  compared  to  the  normalized  AMO  spectrum.  As  expected,  this 
graph  does  approach  zero  at  0.89  pm  .  It  is  also  important  to  note  that  the  spectral  re¬ 
sponse  is  closely  correlated  to  the  spectrum,  thus  the  large  power  output. 


It  is  important  that  this  graph  go  to  zero  in  the  correct  place.  This  may  not  affect 
some  cells  under  AMO  because  only  a  relatively  small  amount  of  energy  is  in  the  spec¬ 
trum  above  where  the  response  goes  to  zero,  but  the  spectrum  from  a  1300K  blackbody, 
however,  contains  most  of  its  energy  at  wavelengths  longer  than  the  cell  should  be  able  to 
respond  to.  If  the  spectral  response  does  not  go  to  zero,  the  model  will  allow  the  cell  to 
respond  to  these  higher  wavelengths  and  give  incorrect  results. 

This  GaAs  cell  model  demonstrates  that  ATLAS  is  capable  of  providing  a  very 
accurate  model  of  a  solar  cell. 
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2.  Gallium  Antimonide 

The  next  cell  built  was  a  gallium  antimonide  (GaSb)  cell.  This  material  is  cur¬ 
rently  one  of  the  best  candidates  for  replacing  the  thermocouples  currently  used  in 
RTG’s.  It  was  difficult  to  find  data  to  compare  with  the  results  of  this  model.  Most  pub¬ 
lications  in  this  area  measure  these  cells  under  concentrators  and  skew  the  efficiency  cal¬ 
culations  to  appear  more  favorable.  One  documented  cell  claimed  to  achieve  approxi¬ 
mately  Voc  =  33  mV  and  Isc  =  34.7  mA  under  AMI. 5  from  a  1  cm2  cell  [Ref.  24].  This 
cell  model  seen  in  Figure  30,  measured  under  AMO  produced  the  output  shown  in  Figure 
31. 


0.2058  pm 
P-  3.72E19 


105  pill 
N-2.1E1S 


35  pm 
N-  7E19 
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Figure  3 1 .  GaSb  Solar  Cell  1- V  Characteristics 


In  addition  to  testing  this  cell  under  AMO,  the  incident  spectrum  was  altered  to 
model  the  spectrum  from  a  concentrator  at  68  and  70  suns  concentration  [Ref .  25].  The 
results  from  these  runs  very  closely  correlated  with  published  data  on  this  cell. 

From  Equation  2.2  the  spectral  response  of  a  GaSb  cell  should  approach  zero  at  a 
wavelength  of  1.7  pm  .  Figure  32  shows  the  spectral  response  of  a  GaSb  cell  compared 
to  the  normalized  AMO  spectrum.  It  can  be  seen  from  the  figure  that  this  response  peaks 
at  X  -  1 .65  pm  ,  and,  as  expected,  it  approaches  zero  at  the  correct  wavelength. 
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3.  Indium  Gallium  Arsenide 

The  final  cell  built  was  an  indium  gallium  arsenide.  The  bandgap  of  this  material 
is  related  to  the  molar  concentration  of  gallium  in  the  compound.  By  varying  this  con¬ 
centration,  the  bandgap  can  be  varied  from  0.35  eV  to  1.42  eV  [Ref.  26].  At  an  In/Ga 
concentration  of  0.72/0.28  InGaAs  has  an  energy  bandgap  of  0.55  eV  [Ref.  14].  It  is  this 
cell  that  was  modeled  here. 

Unlike  GaSb,  there  was  a  great  deal  of  published  data  on  InGaAs  cells  under 
AMO  and  other  spectra  as  well.  Unfortunately,  because  research  on  this  material  is  in  the 
early  stages,  there  is  still  a  great  deal  of  inconsistency  in  published  results.  The  recorded 
short  circuit  current  data  varied  greatly  from  30mA/cm2  to  60mA/cm2  and  recorded 
open- circuit  voltage  data  varied  from  0.2  V  to  0.4  V  ,  so  it  was  difficult  to  determine  the 
accuracy  of  this  model  [Ref.  26].  The  cell  in  Figure  33  was  built  and  found  to  give  re¬ 
sults  that  fell  within  this  range.  The  actual  output  of  this  model  can  be  seen  in  Figure  34. 
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The  spectral  response  of  this  InGaAs  cell  should  approach  zero  at  a  wavelength  of 
2.25  pm  according  to  Equation  2.2.  Figure  35  shows  the  computed  spectral  response  of 

this  cell.  The  spectral  response  peaks  at  X  =  1.55  pm  and  approaches  zero  at  the  correct 
wavelength. 
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Spectral  Response  vs.  AMO 


Figure  35.  Spectral  Response  of  InGaAs 


All  of  the  cells  tested  provided  an  accurate  model  under  the  AMO  spectrum  and 
produced  the  expected  spectral  responses. 


B.  THERMOPHOTO V OLTAIC  CELLS 


The  previous  chapter  established  the  validity  of  the  cell  models.  The  next  step 
was  to  take  these  models  and  optimize  them  for  the  spectrum  of  a  1300K  blackbody.  In 
order  to  do  this,  the  Matlab  programs  previously  discussed  were  used  to  run  iteration  rou¬ 
tines  to  find  the  optimum  thicknesses  and  doping  levels  for  each  layer  of  a  cell. 

1.  Gallium  Antimonide 

Figure  36  below  shows  the  optimum  GaSb  cell  found. 
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Figure  36.  GaSb  TPV  Cell 


Figure  37  shows  the  output  from  this  cell. 


The  photogeneration  rate  was  also  measured  from  this  cell.  Figure  38  shows  the 
photogeneration  rate  at  incrementally  increasing  wavelengths  from  X  =  0.2  pm  -  2.0  pm  . 
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In  this  picture,  brighter  colors  represent  higher  photogeneration  rates.  It  can  be 
seen  that  longer  wavelengths  of  light  are  absorbed  deeper  in  the  cell  and  cause  more  pho- 
togeneration  to  occur.  At  very  short  and  very  long  wavelengths  of  light,  however,  there 
is  almost  no  photogeneration. 


2.  Indium  Gallium  Arsenide 

Figure  39  shows  the  optimum  InGaAs  cell  found. 
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Figure  39.  InGaAs  TPV  Cell 


0.1  jiuii 
N- 1E20 


Figure  40  shows  the  output  of  this  cell  when  the  spectrum  of  1300K  blackbody 
was  applied  to  it. 
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The  results  from  this  cell  were  compared  to  published  data  on  a  similar  InGaAs 
cell  under  a  1225K  blackbody  spectrum. 

The  photogeneration  rate  in  this  cell  was  measured  in  the  same  manner  as  the 
GaSb  cell.  Figure  41  shows  the  photogeneration  for  this  cell  with  incrementally  increas¬ 
ing  wavelengths  of  light  from  X  =  0.2  pm  -  2.3  pm  . 


Again  it  can  be  seen  that  longer  wavelengths  of  light  are  absorbed  deeper  in  the 
cell  and  produce  more  photogeneration  until  the  incident  light  no  longer  has  enough  en¬ 
ergy  to  induce  photogeneration.  From  these  photogeneration  figures  and  the  previous 
spectral  response  graphs,  it  can  be  deduced  that  cells  generally  absorb  light  most  effi¬ 
ciently  that  has  energy  only  slightly  higher  than  the  energy  bandgap  of  the  material. 

This  chapter  compared  simulated  cell  results  under  the  AMO  spectrum  to  meas¬ 
ured  data  in  order  to  verify  the  models.  The  cells  were  then  optimized  for  the  spectrum 
from  a  1300K  blackbody  and  the  results  of  these  cells  were  presented.  The  next  chapter 
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discusses  conclusions  reached  from  this  thesis  and  recommendations  for  future  work  in 
this  area. 
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vn.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  developed  and  demonstrated  reliable  models  of  GaSb  and  InGaAs 
thermophotovoltaic  cells  for  both  the  AMO  spectrum  and  the  spectrum  of  1300K  black- 
body.  In  addition  to  the  Silvaco  models,  several  programs  were  written  to  enable  Matlab 
to  communicate  with  Silvaco  for  mnning  simulations  and  analyzing  data.  With  the  de¬ 
velopment  of  models  for  other  cell  materials  and  possible  multijunction  configurations, 
Silvaco ’s  ATLAS  software  could  prove  to  be  an  invaluable  tool  for  TPV  industry  and  for 
the  cell  designers. 

This  thesis  assumed  the  heat  source  to  be  a  1300K  blackbody,  but  with  very  little 
effort  any  other  spectrum  could  be  applied  to  these  cells.  The  Matlab  program  developed 
here  can  produce  the  spectrum  from  any  temperature  blackbody  over  any  range  of  wave¬ 
lengths,  calculating  the  radiance  at  each  wavelength  along  with  the  total  power  under  the 
curve.  The  optimization  routines  would  mn  exactly  the  same  way. 

The  models  built  here  could  easily  be  configured  to  model  other  cell  materials  as 
well.  The  most  difficult  part  of  building  an  accurate  cell  model  is  providing  accurate  ma¬ 
terial  properties.  These  properties  are  difficult  to  find  for  some  of  the  newer,  more  exotic 
ternary  and  quaternary  materials.  Optical  parameters,  one  of  the  most  important  material 
properties  for  building  an  accurate  model,  can  often  only  be  found  for  wavelengths  of 
light  up  to  1  or  2  pm  but,  for  TPV  applications,  this  data  is  necessary  for  much  longer 
wavelengths  of  light.  Silvaco,  while  a  powerful  simulation  tool,  cannot  accurately  solve 
all  the  equations  to  model  a  cell  without  the  proper  variables. 

Once  initial  cell  models  are  built,  however,  they  can  provide  a  cell  designer  in¬ 
stant  feedback  on  the  impact  of  changes  made  in  any  given  design.  With  the  lack  of  ac¬ 
curate  TPV  cell  models,  the  current  design  process  consists  of  building  and  testing  many 
different  cells  in  order  to  find  the  most  efficient  design.  The  method  presented  in  this 
thesis  allows  a  cell  designer  the  option  of  using  any  configuration  of  any  definable  mate¬ 
rials  with  any  conceivable  doping  schemes.  Other  ongoing  research  is  even  adding  the 
capability  of  predicting  the  impact  of  radiation  on  a  cell.  No  other  model  has  shown  the 
versatility  and  ease  of  design  for  TPV  cells  as  the  Silvaco  software  used  here. 
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The  next  step  in  this  research  is  to  model  a  multijunction  TPV  cell.  Recent  re¬ 
search  by  Green  [Ref.  27]  and  Michalopoulos  [Ref.  28]  has  proven  that  Atlas  is  capable 
of  modeling  and  optimizing  multijunction  solar  cells.  A  multijunction  configuration 
would  have  the  potential  to  greatly  improve  the  heat  to  electrical  conversion  efficiency  of 
the  TPV  cell  as  well. 
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APPENDIX  A:  ATLAS  CODE 


This  appendix  contains  the  Silvaco  input  decks  that  were  used  to  build  the  models 
in  this  thesis  for  GaAs,  GaSb,  and  InGaAs. 


A.  GALLIUM  ARSENIDE 

1.  I-V  Curve 


go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (um) 
set  cellwidth=500 
## 

##  Cap/Contact  (%) 
set  contactpercent=8 
## 

##  Standard  #  of  Divisions 
set  divs=10 
## 

##  Contact 

set  contactthickness= . 1 

## 

##  Cap 

set  capthickness= . 3 
set  capdoping=3el 8 
## 

##  Cathode 

set  cathodethickness= . 2 
set  cathodedoping=l . 5el8 
## 

##  Anode 

set  anodethickness=200 
set  anodedoping=7el7 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness= . 5 
set  bsfdoping=8el8 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=100e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsfL-$bsfthickne ss 
set  bsfdivs=$bsfthickness/$divs 


set  anoL=$bsfH 

set  anoH=$anoL-$anodethickness 
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set  anodivs=$anodethickness/$divs 


set  catL=$anoH 

set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/ $divs 

set  capL=$catH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL- $ contact thickness 
set  contdivs=$contactthickness/$divs 
set  contcenter=$cellwidth/2 
set  contRR=$contcenter+$capwidth/ 2 
set  contLL=$contcenter-$capwidth/ 2 

set  lightY=$contH-5 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdiv 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Caps 

region  num=2 
##Vacuums 
region  num=3 
region  num=4 
y . max=$capL 
#Cathode 
region  num=5 
#Anode 

region  num=6 
#BSF 


material=Gold  x . min=$contLL  x . max=$contRR  y.min=$contH 

material=GaAs  x . min=$contLL  x . max=$contRR  y.min=$capH  y.max=$capL 

material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH  y.max=$capL 
material=vacuum  x . min=$contRR  x . max=$cellwidth  y.min=$contH 

material=GaAs  x.min=0  x . max=$cellwidth  y.min=$catH  y.max=$catL 
material=GaAs  x.min=0  x . max=$cellwidth  y.min=$anoH  y.max=$anoL 


region  num=7  material=GaAs  x.min=0  x . max=$cellwidth  y.min=$bsfH  y.max=$bsfL 


##Electrodes 

electrode  name=cathode  material=Gold  x.min=$contLL  x . max=$contRR  y.min=$contH 
y . max=$contL 

electrode  name=anode  material=Gold  x.min=0  x . max=$cellwidth  y.min=$bsfL 
y . max=$bsf L 
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#Doping 

#  Cap 

doping  uniform  region=2  n.type 

#  Cathode 

doping  uniform  region=5  n.type 

#  Anode 

doping  uniform  region=6  p.type 

#  BSF 


doping  uniform  region=7  p.type 


conc=$capdoping 
conc=$ cathodedoping 
conc=$ anodedoping 
conc=$bsf doping 


#Material  properties 
material  COPT=7e-10 

#  Vacuum 

material  material=Vacuum  real . index=3 . 3  imag.index=0 
$  GaAs 

material  material=GaAs  EG300=1.42  PERMITTIVITY=13 . 1  AFFINITY=4 . 07 
material  material=GaAs  NC300=4.7el7  NV300=7el8 
material  material=GaAs  index . file=GaAs . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 

#  Models 

models  region=3  CONMOB 
models  region=4  CONMOB 
models  region=8  CONMOB 
models  region=9  CONMOB 
models  region=10  CONMOB 
models  OPTR  print 

#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . file=AM0silv. spec  wavel . start=0 . 21  wavel.end=4  wavel.num=50 


###  SOLVING  I-V  CURVE 

set  suns=l 
set  I=$suns* . 95 
set  a=$I* . 9 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=$a 

method  newton  maxtraps=10  itlimit=100 
solve  bl=$I 

extract  name="isc"  max (curve (v. "anode"  ,  i . "cathode" ) ) 

set  isc=$isc*$3D 
set  il=$isc/10 
set  i2=$il+$isc/10 
set  i3=$i2+$isc/10 
set  i4=$i3+$isc/10 
set  i5=$i4+$isc/10 
set  i6=$i5+$isc/20 
set  i7=$i 6+$isc/20 
set  i8=$i7+$isc/20 
set  i9=$i8+$isc/20 
set  il0=$i9+$isc/20 
set  ill=$il0+$isc/40 
set  il2=$il l+$isc/4 0 
set  il3=$il2+$isc/4 0 
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set  il4=$il3+$isc/40 
set  il5=$il4+$isc/40 
set  il6=$il5+$isc/80 
set  il7=$il6+$isc/80 
set  il8=$il7+$isc/80 
set  il9=$il 8+$isc/80 
set  i20=$il 9+$isc/80 
set  i2 l=$i20+$isc/80 
set  i22=$i2 l+$isc/80 
set  i23=$i22+$isc/80 
set  i24=$i23+$isc/80 
set  i25=$i24+$isc/80 

contact  name=anode  current 

method  newton  maxtraps=4  itlimit=100 

log  outfile=gaas . log 
struct  outf ile=gaas . str 

solve  ianode=0  bl=$I 
solve  ianode=-$il  bl=$I 
solve  ianode=-$i2  bl=$I 
solve  ianode=-$i3  bl=$I 
solve  ianode=-$i4  bl=$I 
solve  ianode=-$i5  bl=$I 
solve  ianode=-$i6  bl=$I 
solve  ianode=-$i7  bl=$I 
solve  ianode=-$i8  bl=$I 
solve  ianode=-$i9  bl=$I 
solve  ianode=-$ilO  bl=$I 
solve  ianode=-$ill  bl=$I 
solve  ianode=-$il2  bl=$I 
solve  ianode=-$il3  bl=$I 
solve  ianode=-$il4  bl=$I 
solve  ianode=-$il5  bl=$I 
solve  ianode=-$il6  bl=$I 
solve  ianode=-$il7  bl=$I 
solve  ianode=-$il8  bl=$I 
solve  ianode=-$il9  bl=$I 
solve  ianode=-$i20  bl=$I 
solve  ianode=-$i21  bl=$I 
solve  ianode=-$i22  bl=$I 
solve  ianode=-$i23  bl=$I 
solve  ianode=-$i24  bl=$I 
solve  ianode=-$i25  bl=$I 

log  off 


tonyplot  gaas.log  -set  onecon.set 
tonyplot  gaas.str 


2.  Frequency  Response 

go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (urn) 


74 


set  cellwidth=500 

## 

##  Cap/Contact  (%) 
set  contactpercent=8 
## 

##  Standard  #  of  Divisions 
set  divs=10 
## 

##  Contact 

set  contactthickness= . 1 

## 

##  Cap 

set  capthickness= . 3 
set  capdoping=3el8 
## 

##  Cathode 

set  cathodethickness= . 2 
set  cathodedoping=l . 5el8 
## 

##  Anode 

set  anodethickness=200 
set  anodedoping=7el7 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness= . 5 
set  bsfdoping=8el8 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=100e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsfL-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 

set  anoL=$bsfH 

set  anoH=$anoL-$anodethickness 
set  anodivs=$anodethickness/ $divs 

set  catL=$anoH 

set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/ $divs 

set  capL=$catH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL- $ contact thickness 
set  contdivs=$contactthickness/ $divs 
set  contcenter=$cellwidth/2 
set  contRR=$contcenter+$capwidth/ 2 
set  contLL=$contcenter-$capwidth/ 2 

set  lightY=$contH-5 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
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x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdiv 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Caps 

region  num=2 
##Vacuums 
region  num=3 
region  num=4 
y . max=$capL 
tCathode 
region  num=5 
#Anode 

region  num=6 
#BSF 

region  num=7 


material=Gold  x . min=$contLL  x . max=$contRR  y.min=$contH 

material=GaAs  x.min=$contLL  x . max=$contRR  y.min=$capH  y.max=$capL 

material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH  y.max=$capL 
material=vacuum  x . min=$contRR  x .max=$cellwidth  y.min=$contH 

material=GaAs  x.min=0  x . max=$cellwidth  y.min=$catH  y.max=$catL 
material=GaAs  x.min=0  x . max=$cellwidth  y.min=$anoH  y.max=$anoL 
material=GaAs  x.min=0  x . max=$cellwidth  y.min=$bsfH  y.max=$bsfL 


##Electrodes 

electrode  name=cathode  material=Gold  x . min=$contLL  x.max=$contRR  y.min=$contH 
y . max=$contL 

electrode  name=anode  material=Gold  x.min=0  x . max=$cellwidth  y.min=$bsfL 
y ,max=$bsfL 


#Doping 

#  Cap 

doping  uniform  region=2  n.type 

#  Cathode 

doping  uniform  region=5  n.type 

#  Anode 

doping  uniform  region=6  p.type 

#  BSF 

doping  uniform  region=7  p.type 


conc=$capdoping 
conc=$ cathodedoping 
conc=$anodedoping 
conc=$bsf doping 


#Material  properties 
material  COPT=7e-10 

#  Vacuum 

material  material=Vacuum  real . index=3 . 3  imag.index=0 
$  GaAs 

material  material=GaAs  EG300=1.42  PERMITTIVITY=13 . 1  AFFINITY=4 . 07 
material  material=GaAs  NC300=4.7el7  NV300=7el8 
material  material=GaAs  index . file=GaAs . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 

#  Models 

models  region=3  CONMOB 
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models  region=4  CONMOB 
models  region=8  CONMOB 
models  region=9  CONMOB 
models  region=10  CONMOB 
models  OPTR  print 

#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . file=AMOsilv . spec  wavel . start=0 . 21  wavel.end=4  wavel.num=50 

###  SOLVING  I-V  CURVE 

set  suns=l 
set  I=$suns* . 95 
set  a=$I* . 9 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=$a 

method  newton  maxtraps=10  itlimit=100 
solve  bl=$I 


extract  name="isc"  max (curve (v. "anode" ,  i . "cathode" ) ) 


output  con. band  val.band 
solve  init 

log  outfile=freq. log 


solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0.1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 


lambda=0 . 22 
lambda=0 .225 
lambda=0 . 23 
lambda=0 .235 
lambda=0 . 24 
lambda=0 .245 
lambda=0 . 25 
lambda=0 .255 
lambda=0 . 2  6 
lambda=0 .265 
lambda=0 . 27 
lambda=0 .275 
lambda=0 . 28 
lambda=0 .285 
lambda=0 . 2  9 
lambda=0 .295 
lambda=0 . 30 
lambda=0 .305 
lambda=0 . 31 
lambda=0 . 315 
lambda=0 . 32 
lambda=0 .325 
lambda=0 . 33 
lambda=0 .335 
lambda=0 . 34 
lambda=0 . 345 
lambda=0 . 35 
lambda=0 .355 
lambda=0 . 3  6 
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solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 


lambda=0 . 365 
lambda=0 . 37 
lambda=0 . 375 
lambda=0 . 38 
lambda=0 . 385 
lambda=0 . 3  9 
lambda=0 . 395 
lambda=0 . 4  0 
lambda=0 .405 
lambda=0 . 4 1 
lambda=0 .415 
lambda=0 . 42 
lambda=0 .425 
lambda=0 . 43 
lambda=0 . 435 
lambda=0 . 4  4 
lambda=0 .445 
lambda=0 . 4  5 
lambda=0 . 455 
lambda=0 . 4  6 
lambda=0 . 465 
lambda=0 . 4  7 
lambda=0 .475 
lambda=0 . 4  8 
lambda=0 .485 
lambda=0 . 4  9 
lambda=0 .495 
lambda=0 . 50 
lambda=0 .505 
lambda=0 . 51 
lambda=0 . 515 
lambda=0 . 52 
lambda=0 .525 
lambda=0 . 53 
lambda=0 .535 
lambda=0 . 54 
lambda=0 .545 
lambda=0 . 55 
lambda=0 .555 
lambda=0 . 5  6 
lambda=0 . 565 
lambda=0 . 57 
lambda=0 . 575 
lambda=0 . 58 
lambda=0 . 585 
lambda=0 . 5  9 
lambda=0 . 595 
lambda=0 . 60 
lambda=0 .605 
lambda=0 . 61 
lambda=0 . 615 
lambda=0 . 62 
lambda=0 .625 
lambda=0 . 63 
lambda=0 .635 
lambda=0 . 64 
lambda=0 . 645 
lambda=0 . 65 
lambda=0 .655 
lambda=0 . 66 
lambda=0 . 665 
lambda=0 . 67 
lambda=0 . 675 
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solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 
solve  bl=0 . 1 
solve  bl=0.1 


lambda=0 . 68 
lambda=0 . 685 
lambda=0 . 69 
lambda=0 . 695 
lambda=0 . 7  0 
lambda=0 .705 
lambda=0 . 71 
lambda=0 .715 
lambda=0 . 72 
lambda=0 .725 
lambda=0 . 73 
lambda=0 .735 
lambda=0 . 7  4 
lambda=0 .745 
lambda=0 . 7  5 
lambda=0 .755 
lambda=0 . 7  6 
lambda=0 . 7  65 
lambda=0 . 77 
lambda=0 .775 
lambda=0 . 7  8 
lambda=0 .785 
lambda=0 . 7  9 
lambda=0 .795 
lambda=0 . 80 
lambda=0 . 805 
lambda=0 . 81 
lambda=0 .815 
lambda=0 . 82 
lambda=0 .825 
lambda=0 . 83 
lambda=0 . 835 
lambda=0 . 84 
lambda=0 .845 
lambda=0 . 85 
lambda=0 . 855 
lambda=0 . 8  6 
lambda=0 . 865 
lambda=0 . 87 
lambda=0 .875 
lambda=0 . 88 
lambda=0 .885 
lambda=0 . 8  9 
lambda=0 . 895 
lambda=0 . 90 
lambda=0 .905 
lambda=0 . 91 
lambda=0 . 915 
lambda=0 . 92 
lambda=0 . 925 
lambda=0 . 93 
lambda=0 .935 
lambda=0 . 94 
lambda=0 . 945 
lambda=0 . 95 
lambda=0 .955 
lambda=0 . 96 
lambda=0 . 965 
lambda=0 . 97 
lambda=0 . 975 
lambda=0 . 98 
lambda=0 . 985 
lambda=0 . 99 
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solve  bl=0 . 1  lambda=0.995 
solve  bl=0.1  lambda=1.00 
log  off 


B.  GALLIUM  ANTIMONIDE 

1.  I-V  Curve 

go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (um) 
set  cellwidth=500 
## 

##  Cap/Contact  (%) 
set  contactpercent=8 
## 

##  Standard  #  of  Divisions 
set  divs=20 
## 

##  Contact 

set  contactthickness= . 1 

## 

##  Cap  (there  is  no  cap  in  this  model) 
set  capthickness= . 1 
set  capdoping=3el8 
## 

##  Cathode (N-type ) 

set  cathodethickness=10 

set  cathodedoping=l . 000000e+020 

## 

##  Anode (P-type )  (.2-.5um) 

set  anodethickness= . 37 
set  anodedoping=9 . 8el9 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness=15 . 000000 
set  bsfdoping=l . 000000e+020 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=100e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsf L-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 

set  catL=$bsfH 

set  catH=$catL-$cathodethickness 


80 


set  catdivs=$cathodethickness/ $divs 
set  anoL=$catH 

set  anoH=$anoL-$anodethickness 
set  anodivs=$anodethickness/ $divs 

set  capL=$anoH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL-$ contact thickness 
set  contdivs=$contactthickness/ $divs 
set  contcenter=$cellwidth/2 
set  contRR=$contcenter+$capwidth/ 2 
set  contLL=$contcenter-$capwidth/ 2 

#Sets  light  source  200um  above  surface  (try  different  heights) 
set  lightY=$contH-200 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdiv 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Cap 

region  num=2 
y . max=$capL 
##Vacuums 
region  num=3 
y . max=$capL 
region  num=4 
y . min=$contH 
#Anode 

region  num=5 
y . max=$anoL 
#Cathode 
region  num=6 
y . max=$catL 
#BSF 


material=Gold  x . min=$contLL  x.max=$contRR  y.min=$contH 


material=GaSb  x . min=$contLL  x . max=$contRR  y.min=$capH 


material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH 

material=vacuum  x . min=$contRR  x . max=$cellwidth 
y . max=$capL 

material=GaSb  x.min=0  x . max=$cellwidth  y.min=$anoH 

material=GaSb  x.min=0  x . max=$cellwidth  y.min=$catH 
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region  num=7  material=GaSb  x.min=0  x .max=$cellwidth  y.min=$bsfH 
y . max=$bsf L 

##Electrodes 

electrode  name=anode  material=Gold  x . min=$contLL  x.max=$contRR 
y.min=$contH  y.max=$contL 

electrode  name=cathode  material=Gold  x.min=0  x . max=$cellwidth 
y.min=$bsfL  y.max=$bsfL 


#Doping  (In  this  model 
eral ) 

#  Cap 

doping  uniform  region=2 

#  Anode 

doping  uniform  region=5 

#  Cathode 

doping  uniform  region=6 

#  BSF 


doping  uniform  region=7 


there  are  no  specified  doping 

p.type  conc=$capdoping 
p.type  conc=$anodedoping 
n.type  conc=$cathodedoping 
n.type  conc=$bsfdoping 


levels , 


try  sev- 


#Material  properties 

#  Vacuum 

material  material=Vacuum  real . index=3 . 8  imag.index=0 

#  GaSb 

material  material=GaSb  COPT=5e-10 

material  material=GaSb  EG300=.726  PERMITTIVITY=14 . 4  AFFINITY=4 . 06 
material  material=GaSb  NC300=2.1el7  NV300=1.8el9 
material  material=GaSb  index . file=GaSb . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 

#  Models 

models  OPTR  print 
#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . file=1300BB . spec  wavel . start=0 . 207  wavel.end=5  wavel . num=1000 


###  SOLVING  I-V  CURVE 

set  suns=l 
set  I=$suns* . 95 
set  a=$ I  *  .  9 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=$a 

method  newton  maxtraps=10  itlimit=100 
solve  bl=$I 

extract  name="isc"  max (curve (v anode " ,  i cathode ") ) 


set  isc=$isc*$3D 
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set  il=$isc/10 
set  i2=$il+$isc/10 
set  i3=$i2+$isc/10 
set  i4=$i3+$isc/10 
set  i5=$i4+$isc/10 
set  i6=$i5+$isc/20 
set  i7=$i6+$isc/20 
set  i8=$i7+$isc/20 
set  i9=$i8+$isc/20 
set  il0=$i9+$isc/20 
set  ill=$il0+$isc/40 
set  il2=$ill+$isc/40 
set  il3=$il2+$isc/40 
set  il4=$il3+$isc/40 
set  il5=$il4+$isc/40 
set  il6=$il5+$isc/80 
set  il7=$il6+$isc/80 
set  il8=$il7+$isc/80 
set  il9=$il8+$isc/80 
set  i20=$il9+$isc/80 
set  i21=$i20+$isc/80 
set  i22=$i21+$isc/80 
set  i23=$i22+$isc/80 
set  i24=$i23+$isc/80 
set  i25=$i24+$isc/80 

contact  name=anode  current 

method  newton  maxtraps=4  itlimit=100 

log  outf ile=gasb . log 
struct  outf ile=gasb . str 

solve  ianode=0  bl=$I 
solve  ianode=-$il  bl=$I 
solve  ianode=-$i2  bl=$I 
solve  ianode=-$i3  bl=$I 
solve  ianode=-$i4  bl=$I 
solve  ianode=-$i5  bl=$I 
solve  ianode=-$i6  bl=$I 
solve  ianode=-$i7  bl=$I 
solve  ianode=-$i8  bl=$I 
solve  ianode=-$i9  bl=$I 
solve  ianode=-$ilO  bl=$I 
solve  ianode=-$ill  bl=$I 
solve  ianode=-$il2  bl=$I 
solve  ianode=-$il3  bl=$I 
solve  ianode=-$il4  bl=$I 
solve  ianode=-$il5  bl=$I 
solve  ianode=-$il6  bl=$I 
solve  ianode=-$il7  bl=$I 
solve  ianode=-$il8  bl=$I 
solve  ianode=-$il9  bl=$I 
solve  ianode=-$i20  bl=$I 
solve  ianode=-$i21  bl=$I 
solve  ianode=-$i22  bl=$I 
solve  ianode=-$i23  bl=$I 
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solve  ianode=-$i24  bl=$I 
solve  ianode=-$i25  bl=$I 


log  off 

tonyplot  gasb.log  -set  onecon.set 
tonyplot  gasb.str 

2.  Frequency  Response 

go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (um) 
set  cellwidth=500 
## 

##  Cap/Contact  (%) 
set  contactpercent=8 
## 

##  Standard  #  of  Divisions 
set  divs=20 
## 

##  Contact 

set  contactthickness= . 1 

## 

##  Cap  (there  is  no  cap  in  this  model) 
set  capthickness= . 1 
set  capdoping=3el 8 
## 

##  Cathode (N-type ) 

set  cathodethickness=10 

set  cathodedoping=l . 000000e+020 

## 

##  Anode (P-type )  (.2-.5um) 

set  anodethickness= . 37 
set  anodedoping=9 . 8el9 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness=15 . 000000 
set  bsfdoping=l . 000000e+020 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=50e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsfL-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 

set  catL=$bsfH 

set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/ $divs 

set  anoL=$catH 

set  anoH=$anoL-$anodethickness 
set  anodivs=$anodethickness/$divs 
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set  capL=$anoH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL- $ contact thickness 
set  contdivs=$contactthickness/$divs 
set  contcenter=$cellwidth/2 
set  contRR=$contcenter+$capwidth/ 2 
set  contLL=$contcenter-$capwidth/ 2 

#Sets  light  source  200um  above  surface  (try  different  heights) 
set  lightY=$contH-200 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdi v 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Cap 

region  num=2 
##Vacuums 
region  num=3 
region  num=4 
y . max=$capL 
#Anode 

region  num=5 
tCathode 
region  num=6 
#BSF 

region  num=7 


material=Gold  x . min=$contLL  x ,max=$contRR  y.min=$contH 

material=GaSb  x.min=$contLL  x . max=$contRR  y.min=$capH  y.max=$capL 

material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH  y.max=$capL 
material=vacuum  x . min=$contRR  x ,max=$cellwidth  y.min=$contH 

material=GaSb  x.min=0  x . max=$cellwidth  y.min=$anoH  y.max=$anoL 
material=GaSb  x.min=0  x . max=$cellwidth  y.min=$catH  y.max=$catL 
material=GaSb  x.min=0  x . max=$cellwidth  y.min=$bsfH  y.max=$bsfL 


##Electrodes 

electrode  name=anode  material=Gold  x . min=$contLL  x.max=$contRR  y.min=$contH 
y . max=$contL 

electrode  name=cathode  material=Gold  x.min=0  x . max=$cellwidth  y.min=$bsfL 
y .max=$bsfL 


#Doping  (In  this  model 

#  Cap 

doping  uniform  region=2 

#  Anode 

doping  uniform  region=5 

#  Cathode 

doping  uniform  region=6 


there  are  no  specified  doping 
p.type  conc=$capdoping 
p.type  conc=$anodedoping 
n.type  conc=$cathodedoping 


levels , 


try  several) 
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#  BSF 

doping  uniform  region=7  n.type  conc=$bsfdoping 


material 

#  Vacuum 
material 

#  GaSb 
material 
material 
material 

#  Gold 


C0PT=7e- 10 

material=Vacuum  real . index=3 . 3  imag.index=0 

material=GaSb  EG300=.72  PERMITTIVITY=15 . 7  AFFINITY=4 . 06 
material=GaSb  NC300=2 . 07 84 6el 7  NV300=1 . 818 65el7 
material=GaSb  index . f ile=GaAs . opt 


material  material=Gold  real . index=l . 2  imag . index=l . 8 
#  Models 

models  OPTR  print 


#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 
power . file=1300BB . spec  wavel . start=0 . 1  wavel.end=7  wavel . num=100 


set  1= . 1 


solve  init 
solve  bl=$I 

extract  name="isc"  max (curve (v. "anode"  ,  i . "cathode" ) ) 

output  con. band  val.band 
solve  init 

log  outfile=freq. log 
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C.  INDIUM  GALLIUM  ARSENIDE 


1.  I-V  Curve 

go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (um) 
set  cellwidth=500 
## 

##  Cap/Contact  (%) 
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set  contactpercent=8 

## 

##  Standard  #  of  Divisions 
set  divs=20 
## 

##  Contact 

set  contactthickness= .  1 

## 

## 

##Molar  Concentration  of  In ( 1-x) Ga (x) As (x) 
set  x=.28 

set  perm=15 . 1-2 . 87*$x+ . 67*$xA2 
set  af f in=4 . 9- . 83*$x 
set  n=3 . 5 1- . 2 1  * $x 
## 

##  Cap  (there  is  no  cap  in  this  model) 
set  capthickness= . 1 
set  capdoping=3el8 
## 

##  Cathode (N-type ) 
set  cathodethickness=17 
set  cathodedoping=le20 
## 

##  Anode (P-type )  (.2-.5um) 

set  anodethickness=2 . 22 
set  anodedoping=le20 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness= . 1 
set  bsfdoping=le20 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=100e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsf L-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 

set  catL=$bsfH 

set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/ $divs 

set  anoL=$catH 

set  anoH=$anoL-$anodethickness 
set  anodivs=$anodethickness/ $divs 

set  capL=$anoH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL-$ contact thickness 
set  contdivs=$contactthickness/ $divs 
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set  contcenter=$cellwidth/2 

set  contRR=$contcenter+$capwidth/ 2 

set  contLL=$contcenter-$capwidth/ 2 

#Sets  light  source  200um  above  surface  (try  different  heights) 
set  lightY=$contH-200 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdiv 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Cap 

region  num=2 
y . max=$capL 
##Vacuums 
region  num=3 
y . max=$capL 
region  num=4 
y . min=$contH 
#Anode 

region  num=5 
y . max=$anoL 
#Cathode 
region  num=6 
y . max=$catL 
#BSF 


material=Gold  x . min=$contLL  x.max=$contRR  y.min=$contH 


material=InGaAs  x . min=$contLL  x.max=$contRR  y.min=$capH 

material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH 

material=vacuum  x . min=$contRR  x . max=$cellwidth 
y . max=$capL 

material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$anoH 
material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$catH 


region  num=7  material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$bsfH 
y . max=$bsf L 


##Electrodes 

electrode  name=anode  material=Gold  x . min=$contLL  x.max=$contRR 
y.min=$contH  y.max=$contL 

electrode  name=cathode  material=Gold  x.min=0  x . max=$cellwidth 
y.min=$bsfL  y.max=$bsfL 


#Doping  (In  this  model  there  are  no  specified  doping  levels,  try  sev¬ 
eral  ) 
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#  Cap 

doping  uniform 

#  Anode 

doping  uniform 

#  Cathode 
doping  uniform 

#  BSF 

doping  uniform 


region=2  p.type 
region=5  p.type 
region=6  n.type 
region=7  n.type 


conc=$ capdoping 
conc=$anodedoping 
conc=$ cathodedoping 
conc=$bsf doping 


#Material  properties 

#  Vacuum 

material  material=Vacuum  real . index=3 . 4  imag.index=0 

#  InGaSb 

material  material=InGaAs  COPT=7e-10 

material  material=InGaAs  EG300=.55  PERMI TT IVI TY=$perm  AFFINITY=$af f in 
material  material=InGaAs  NC300=1 . 15el7  NV300=8 . 12el8 
material  material=InGaAs  index . file=InGaAs . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 

#  Models 

models  OPTR  print 
#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 
power . file=1300bb . spec  wavel . start=0 . 2 1  wavel.end=7  wavel.num=50 


###  SOLVING  I-V  CURVE 

set  suns=l 
set  1= . 95*$suns 
set  a=$ I  * . 75 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=$a 

method  newton  maxtraps=10  itlimit=100 
solve  bl=$I 

extract  name="isc"  max (curve (v ." anode " ,  i ." cathode ") ) 

set  isc=$isc*$3D 
set  il=$isc/10 
set  i2=$il+$isc/10 
set  i3=$i2+$isc/10 
set  i4=$i3+$isc/10 
set  i5=$i4+$isc/10 
set  i6=$i5+$isc/20 
set  i7=$i6+$isc/20 
set  i8=$i7+$isc/20 
set  i9=$i8+$isc/20 
set  il0=$i9+$isc/20 
set  ill=$il0+$isc/40 
set  il2=$ill+$isc/40 
set  il3=$il2+$isc/40 
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set  il4=$il3+$isc/40 
set  il5=$il4+$isc/40 
set  il6=$il5+$isc/80 
set  il7=$il6+$isc/80 
set  il8=$il7+$isc/80 
set  il9=$il8+$isc/80 
set  i20=$il9+$isc/80 
set  i21=$i20+$isc/80 
set  i22=$ i2 l+$ isc/ 8 0 
set  i23=$i22+$isc/80 
set  i24=$i23+$isc/80 
set  i25=$i24+$isc/80 

contact  name=anode  current 

method  newton  maxtraps=10  itlimit=100 

log  outf ile=ingaas . log 
struct  outf ile=ingaas . str 

solve  ianode=0  bl=$I 
solve  ianode=-$il  bl=$I 
solve  ianode=-$i2  bl=$I 
solve  ianode=-$i3  bl=$I 
solve  ianode=-$i4  bl=$I 
solve  ianode=-$i5  bl=$I 
solve  ianode=-$i6  bl=$I 
solve  ianode=-$i7  bl=$I 
solve  ianode=-$i8  bl=$I 
solve  ianode=-$i9  bl=$I 
solve  ianode=-$ilO  bl=$I 
solve  ianode=-$ill  bl=$I 
solve  ianode=-$il2  bl=$I 
solve  ianode=-$il3  bl=$I 
solve  ianode=-$il4  bl=$I 
solve  ianode=-$il5  bl=$I 
solve  ianode=-$il6  bl=$I 
solve  ianode=-$il7  bl=$I 
solve  ianode=-$il8  bl=$I 
solve  ianode=-$il9  bl=$I 
solve  ianode=-$i20  bl=$I 
solve  ianode=-$i21  bl=$I 
solve  ianode=-$i22  bl=$I 
solve  ianode=-$i23  bl=$I 
solve  ianode=-$i24  bl=$I 
solve  ianode=-$i25  bl=$I 

log  off 

tonyplot  ingaas.log  -set  onecon.set 
tonyplot  ingaas.str 
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2.  Frequency  Response 

go  atlas 

############################################# 
##  Parameters 
## 

##  Cell  Width  (um) 
set  cellwidth=500 
## 

##  Cap/Contact  (%) 
set  contactpercent=8 
## 

##  Standard  #  of  Divisions 
set  divs=20 
## 

##  Contact 

set  contactthickness= . 1 

## 

## 

##Molar  Concentration  of  In (1-x) Ga (x) As (x) 
set  x=.28 

set  perm=15 . 1-2 . 87*$x+ . 67*$xA2 
set  af f in=4 . 9- . 83*$x 
set  n=3 . 5 1- .  2 1  * $x 

##  Cap  (there  is  no  cap  in  this  model) 
set  capthickness= . 1 
set  capdoping=3el8 
## 

##  Cathode (N-type ) 
set  cathodethickness=17 . 000000 
set  cathodedoping=le20 
## 

##  Anode (P-type )  (.2-.5um) 

set  anodethickness=2 . 22 
set  anodedoping=l+20 
## 

##  Back  Surface  Field  (BSF) 
set  bsfthickness= . 1 
set  bsfdoping=le20 
## 

############################################# 

set  cellwdiv=$cellwidth/$divs 
set  3D=100e6  /$cellwidth 

set  capwidth= . 01*$contactpercent*$cellwidth 
set  capwdiv=$capwidth/$divs 

set  bsfL=0 

set  bsfH=$bsf L-$bsfthickness 
set  bsfdivs=$bsfthickness/$divs 

set  catL=$bsfH 

set  catH=$catL-$cathodethickness 
set  catdivs=$cathodethickness/ $divs 
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set  anoL=$catH 

set  anoH=$anoL-$anodethickness 
set  anodivs=$anodethickness/ $divs 

set  capL=$anoH 

set  capH=$capL-$capthickness 
set  capdivs=$capthickness/$divs 

set  contL=$capH 

set  contH=$contL-$ contact thickness 
set  contdivs=$contactthickness/ $divs 
set  contcenter=$cellwidth/2 
set  contRR=$contcenter+$capwidth/ 2 
set  contLL=$contcenter-$capwidth/ 2 

#Sets  light  source  200um  above  surface  (try  different  heights) 
set  lightY=$contH-200 


mesh  width=$3D 
#X-Mesh 

x.mesh  loc=0  spac=$cellwdiv 
x.mesh  loc=$contLL  spac=$capwdiv 
x.mesh  loc=$contRR  spac=$capwdiv 
x.mesh  loc=$cellwidth  spac=$cellwdiv 


#Y-Mesh 

y.mesh  loc=$contH  spac=$contdivs 
y.mesh  loc=$capH  spac=$capdivs 
y.mesh  loc=$catH  spac=$catdivs 
y.mesh  loc=$anoH  spac=$anodivs 
y.mesh  loc=$bsfH  spac=$bsfdivs 
y.mesh  loc=$bsfL  spac=$bsfdivs 


#Regions 
##Contacts 
region  num=l 
y . max=$contL 
##Cap 

region  num=2 
y . max=$capL 
##Vacuums 
region  num=3 
y . max=$capL 
region  num=4 
y . min=$contH 
#Anode 

region  num=5 
y . max=$anoL 
#Cathode 
region  num=6 
y . max=$catL 
#BSF 

region  num=7 
y . max=$bsf L 


material=Gold  x . min=$contLL  x.max=$contRR  y.min=$contH 
material=InGaAs  x . min=$contLL  x . max=$contRR  y.min=$capH 

material=Vacuum  x.min=0  x . max=$contLL  y.min=$contH 

material=vacuum  x . min=$contRR  x . max=$cellwidth 
y . max=$capL 

material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$anoH 

material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$catH 

material=InGaAs  x.min=0  x . max=$cellwidth  y.min=$bsfH 
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##Electrodes 

electrode  name=anode  material=Gold  x . min=$contLL  x.max=$contRR 
y.min=$contH  y.max=$contL 

electrode  name=cathode  material=Gold  x.min=0  x . max=$cellwidth 
y.min=$bsfL  y.max=$bsfL 

#Doping  (In  this  model  there  are  no  specified  doping  levels,  try  sev¬ 
eral  ) 

#  Cap 

doping  uniform  region=2  p.type  conc=$capdoping 

#  Anode 

doping  uniform  region=5  p.type  conc=$anodedoping 

#  Cathode 

doping  uniform  region=6  n.type  conc=$cathodedoping 

#  BSF 

doping  uniform  region=7  n.type  conc=$bsfdoping 

#Material  properties 

#  Vacuum 

material  material=Vacuum  real . index=3 . 4  imag.index=0 

#  InGaSb 

material  material=InGaAs  COPT=7e-10 

material  material=InGaAs  EG300=.55  PERMI TT IVI TY=$perm  AFFINITY=$af f in 
material  material=InGaAs  NC300=1 . 15el7  NV300=8 . 12el8 
material  material=InGaAs  index . file=InGaAs . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 

#  Models 

models  OPTR  print 
#Light  beams 

beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . file=1300bb . spec  wavel . start=0 . 2 1  wavel.end=5  wavel . num=100 

###  SOLVING  I-V  CURVE 

set  1= . 1 

solve  init 
solve  bl=$I 

extract  name="isc"  max (curve (v ." anode " ,  i ." cathode ") ) 

output  con. band  val.band 
solve  init 

log  outf ile=f req. log 
solve  bl=0.1  lambda=0.22 
solve  bl=0.1  lambda=0.3 
solve  bl=0.1  lambda=0.4 
solve  bl=0.1  lambda=0.5 
solve  bl=0.1  lambda=0.6 
solve  bl=0.1  lambda=0.7 
solve  bl=0.1  lambda=0.8 
solve  bl=0.1  lambda=.82 
solve  bl=0.1  lambda=.84 
solve  bl=0.1  lambda=.86 
solve  bl=0.1  lambda=.88 
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solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

solve 

bl  =  0 

log  off 

1  lambda=0 . 9 
1  lambda=l . 0 
1  lambda=l .  1 
1  lambda=l . 2 
1  lambda=l .  3 
1  lambda=l . 4 
1  lambda=l .  5 
1  lambda=l . 6 
1  lambda=l .  7 
1  lambda=l . 8 
1  lambda=l .  9 
1  lambda=2 . 0 
1  lambda=2 . 1 
1  lambda=2 . 2 
1  lambda=2 . 3 
1  lambda=2 . 4 
1  lambda=2 . 5 
1  lambda=2 . 6 
1  lambda=2 . 7 
1  lambda=2 . 8 
1  lambda=2 . 9 
1  lambda=3 . 0 
1  lambda=3 . 1 
1  lambda=3 . 2 
1  lambda=3 . 3 
1  lambda=3 . 4 
1  lambda=3 . 5 
1  lambda=3 . 6 
1  lambda=3 . 7 
1  lambda=3 . 8 
1  lambda=3 . 9 
1  lambda=4 . 0 
1  lambda=4 . 1 
1  lambda=4 . 2 
1  lambda=4 . 3 
1  lambda=4 . 4 
1  lambda=4 . 5 
1  lambda=4 . 6 
1  lambda=4 . 7 
1  lambda=4 . 8 
1  lambda=4 . 9 
1  lambda=5 . 0 
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APPENDIX  B:  MATLAB  CODE 


This  section  contains  the  Matlab  functions  that  were  written  and  used  for  this  the¬ 


sis. 


A.  ATLASARUN.M 


%atlasrun . m 

o. 

o 

%This  program  runs  iterations  of  Atlas.  It  uses  ivmaxp.m  to  calculate 
open  %circuit  voltage,  short  circuit  current,  max  power,  fill  factor, 
and  efficiency. 

%The  function  filerw.m  is  called  to  iterate  the  ATLAS  input  deck.  The 
function  %time.m  is  called  to  calculate  the  time  of  each  iteration  in 
hours,  minutes,  and  %seconds. 

o. 

o 

% (c)  2004  by  B.P.  Davenport 

close  all 
clear  all 
clc 

%specifies  the  number  of  iterations  to  perform 
iterations=250; 

%specifies  the  range  of  values  to  iterate  between 
old=l . 8 ; ol=old; 
f inal=5 ; 

step= (final-old) /iterations; 
k=0 ; 

total runt ime=0 ; 

for (a=ol : step : final) 
tic 
k=k+l; 

f idres=f open ( 'atres.txt 1 , ' a ' ) 
f ilerw ( ' ingaas . in ' , old, a) 
delete  done.log 
delete  atlaslog.log 

%Here  atlas  is  called  by  Matlab  to  run  the  specified  input  deck 
! D : \Silvaco\etc\GuiAppStarter . exe  -lib-dir-name  deckbuild  -exe-name 
Deckbld  -run  ingaas . in  -outfile  atlaslog.log 

sprintf (' Executing  run  %u\nStandby  for  results', k) 
xy=-l; 

while (xy==-l ) 

xy=f open ( ' done . log ' ) ; 
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pause  ( 3 ) 


end 

%Here  the  * .  log  file  created  by  ATLAS  is  read  and  open-circuit 
voltage, 

%short  circuit  current,  max  power,  fill  factor  and  efficiency  are 
returned . 

[Voc, I  sc, Pmax, FF, Ef f ] =ivmaxp ( ' ingaas . log ' )  ; 
data(k,  :)  =  [a  Voc  Isc  Pmax  FF  Eff ]  ; 

%The  previously  calculated  data  is  written  to  file, 
fprintf (fidres,  ' %f  %f  %f  %f  %f  %f \n ' , data (k,  : ) )  ; 

%After  each  iteration,  all  calculated  data  is  printed  to  the  screen 
%including  the  time  it  took  to  run  the  iteration.  A  plot  of  effi¬ 
ciency  is 

%updated  after  each  iteration. 

clc 

data 

thisruntime=toc; Tl=time (toe)  ; 

totalruntime=totalruntime+thisruntime; T2=time (totalruntime)  ; 
averageruntime=totalruntime/k; T3=time (averageruntime) ; 
estimaterun- 

time=averageruntime* (iterations  +  1 ) ; T4=time (estimateruntime)  ; 

sprint f (' This  run  took  %s\nTotal  run  time  so  far  is  %s\nAverage  run 
time  is  %s\nEstimated  total  run  time  is  %s ' , Tl, T2, T3, T4) 

old=a; 

f close ( ' all ' ) 
pause ( 3 ) 

plot (data ( : , 1 ) , data ( : , 6 ) ) 
title ( 'Efficiency' ) 

end 


figure ( 1 ) 

plot (data ( : , 1 ) , data ( : , 2 ) ) 
title ( 'Voc ' ) 

figure (2 ) 

plot (data ( : , 1 ) , data ( : , 3 ) ) 
title (' Isc ' ) 

figure  (3) 

plot (data ( : , 1 ) , data  ( : , 4 ) ) 
title ( ' Pmax ' ) 

figure  ( 4 ) 

plot (data ( : , 1 ) , data  ( : , 5 ) ) 
title ( 'FF  '  ) 

figure  (5) 

plot (data ( : , 1 ) , data  ( : , 6 ) ) 
title (' Eff ' ) 
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B 


FILERW.M 


function  filerw (file, old, new) 

%f ilerw (file, old, new) 

%This  program  opens  an  infile,  "file",  and  writes  over  the 
%"old"  anode  thickness  with  the  "new"  anodethickness 
%file  must  be  in  11 

%old  and  new  are  the  old  and  new  values  of  anodethickness 
%This  same  function  was  also  used  to  iterate  other  variables  by  replac¬ 
ing  % 'anodethickness'  with  the  name  of  the  variable  to  be  iterated. 

o. 

o 

% (c)  2004  B.P.  Davenport 

f idr=fopen (file, ' r ' ) ; 
a=f scanf (f idr,  ' %c ' )  ; 
fclose (f idr) ; 
f idw=f open (file, ' w ' ) ; 

a=strrep (a,  sprintf ( ' anodethickness=%f ' , old) , sprintf ( ' anodethickness=%f ' 

,  new) )  ; 

fwrite (fidw, a) ; 
fclose (fidw) ; 


C.  IVMAXP.M 


function  [Voc, Isc, Pmax, FF, Ef f ] =ivmaxp (varargin) 

% [Voc, Isc, Pmax, FF, Ef f ] =ivmaxp (file, type) 

%Finds  Voc,  Isc,  Pmax,  FF,  and  Efficiency  from  an  Atlas  * . log  file 
%file  and  type  must  be  in  '  ' 

%Type  is  the  doping  of  the  top  layer,  n  or  p 

o. 

o 

% (c)  2004  B.P.  Davenport 


if (nargin==l ) 

file=char (varargin  (1) ) ; 
type= ' P  '  ; 
elseif (nargin==2) 

file=char (varargin  (1) ) 
type=char (varargin  (2) ) 

else 

'Type  mus  be  either  "N",  "n",  "P",  or  "p"' 

Voc=NaN ; 

Isc=NaN; 

Pmax=NaN ; 

FF=NaN ; 
return 

end 


if (type== ' N ' | type== ' n ' ) 

[I  V] =textread (file, ' %*s  %*f  %*f  %*f  %*f  %*f  %*f  %f  %f  %*f 

%*f ' ,  ' header lines  '  , 20 )  ; 
elseif (type== ' P ' | type== ' p ' ) 
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[V  I ] =textread ( f ile ,  '%*s  %*f  %*f  %*f  %*f  %f  %*f  %*f  %*f  %*f 

%f ' , ' header lines ' , 20 ) ; 
else 

'Type  mus  be  either  "N",  "n",  "P",  or  "p"' 

Voc=NaN ; 

Isc=NaN; 

Pmax=NaN ; 

FF=NaN ; 
return 

end 

P=I . *V ; 

[a  b] =max  (P) ; 

Pmax=V (b) *1  (b) ; 

Voc=max (V) ; 

I sc=max ( I ) ; 

FF=Pmax/ (Voc*Isc) ; 

Ef f=100*Pmax/ ( . 1353) ; 


x=l . 1 *Voc ; 
y=l . l*Isc; 

plot(V,I,V(b),I(b),  '  rh  '  ) 
xlabel ( 'Voltage  (V) ' ) 
ylabel ( 'Current  (A) ' ) 
title ( ' I-V  Curve') 
axis ( [ 0  x  0  y ] ) 


D.  BLACKBODY.  M 


function  [Lx, Sx] =blackbody (spec) 

% [Lx, Sx] =blackbody (spec) 

%This  program  calculates  the  spectrum  a  5800K  blackbody  which  repre¬ 
sents  AMO  and 

%the  spectrum  from  any  other  temperature  body  specified  by  'T'  .  The 
wavelength  %and  radiance  are  returned  for  the  desired  spectrum.  The 
spectrum  of  the  %blackbody  source  is  written  to  a  * . spec  file  formatted 
for  ATLAS  and  the  %results  are  plotted. 

o. 

o 

% (c)  2004  B.P.  Davenport 


spec=char (spec) ; 

Tsun=5800 ; 

T=1300 ; 

%constants  definition 
h=6 . 6260693e-34; 
c=299792458; 
k=l . 386505e-23 ; 

%scaling  factor  to  account  for  the  distance  from  the  sun 
dilut ion=l / 46200; 
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%the  number  of  evenly  spaced  points  to  calculate 
points=3000 ; 

L=linspace (le-18, lle-6, points) *10A6; 

Ssun=  (2 . *pi . *c . A2 . *h)  ,/((L*10A-6)  .  A5  .  * (exp (h.*c./((L*10A-6)  . *k. *Tsun) ) - 
1) ) *le-10*dilution; 

[xsun,  ysun] =max (Ssun) ; 

S= (2 . *pi . *c. A2 . *h) ./ ( (L*10A-6) .  A5.* (exp(h.*c./ ( (L*10A-6) . *k.*T) )- 
1) ) *le-10; 

[ xot her , y other ] =max (S) ; 

f id=f open ('1300 BB. spec',  ' w ' ) ; 
fprintf(fid, ' %u\n ' , length (S) ) 
for (x=l :points) 

fprintfffid, ' %f  %f\n',L(x),S(x)); 

end 

fclose (fid) 

plot (L, Ssun, ' b ' , L, S, ' g ' ) 

title (' Optical  File') 

xlabel (' Wavelength  (um) ') 

ylabel  (' Power  Density  (W/cmA2/um)  ') 


figure  (2 ) 

plot (L, Ssun, ' b ' ) 

title ('Sun  Optical  File') 

xlabel (' Wavelength  (um) ') 

ylabel (' Power  Density  (W/cmA2/um) ') 

psun=0 ; 
pother=0 ; 

for  (x=l :  (points-1) ) 

psun=psun+ (L (x+1) -L (x) ) *Ssun (x) ; 
pother=pother+ (L (x+1) -L (x) ) *S (x)  ; 

end 

pother 

psun 

pother 

ps  =  sprintf  (' Total  power  is:\n%f  Watts/cmA2\nThe  peak  is:\n%f 
Watts/cmA2\nAt  L=%f  um ', pother, xother, L (yother) ) ; 

ps2=sprintf (' Total  power  is:\n%f  Watts/cmA2\nThe  peak  is:\n%f 
Watts/cmA2\nAt  L=%f  um ', psun, xsun, L (ysun) )  ; 


figure  (3) 
plot (L, S, ' g ' ) 

title ('1300  Blackbody  Radiance') 
xlabel (' Wavelength  (um) ') 
ylabel  (' Power  Density  (W/cmA2/um)  ') 
text ( 1 . 5 , 1 , ps ) 

figure  (2 ) 
text  (2 ,  . 1 , ps2 ) 

Lx=L  '  ; 

if ( strcmp (spec, ' BB ' ) | strcmp ( spec, ' bb ' ) ) 
Sx=S  '  ; 
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end 

if (strcmp (spec,  ' Sun ' )  | strcmp (spec,  ' sun ' ) ) 
Sx=Ssun ' ; 

end 


E.  FREQRESP.M 


function  [L, photo] =f reqresp (varargin) 

o. 

o 

% [L, photo] =f reqresp (varargin) 

%This  program  reads  in  data  a  specified  file  and  returns  and  plots  the 
frequency  %response  from  a  cell. 

o. 

o 

% (c)  2004  B.P.  Davenport 

%If  there  is  only  one  input,  the  cell  is  assumed  to  be  P  on  N. 
if (nargin==l ) 

file=char (varargin  (1) ) ; 
type= ' P  '  ; 
elseif (nargin==2) 

file=char (varargin  (1) ) ; 
type=char (varargin  (2) ) ; 
elseif (nargin==3) 

file=char (varargin  (1) ) ; 
type==char (varargin  (2) ) ; 

else 

'Type  mus  be  either  "N",  "n",  "P",  or  "p"' 

L=NaN; 
photo=NaN ; 
return 

end 

%This  function  can  accept  data  in  two  different  forms  from  ATLAS  or 
from  a  file  %containing  wavelength  radiance  pairs, 
if (type== ' N ' | type== ' n ' ) 

[L  photo]  =textread  ( file,  '%*s  %*f  %*f  %*f  %f  %*f  %*f  %f  %*f  %*f 

%*f '  ,  ' header lines  '  , 20 )  ; 
elseif (type== ' P  '  | type== ' p ' ) 

[L  photo]  =textread  (file,  '%*s  %*f  %*f  %*f  %f  %*f  %*f  %*f  %*f  %*f 

%f ' ,  ' header lines ' , 2  0 )  ; 
elseif (type== ' x ' | type== ' X ' ) 

[L  photo] =t ext re ad (file,  ' %f%f ' ,  'delimiter ' ,  ' \t  '  ) ; 

else 

'Type  mus  be  either  "N",  "n",  "P",  or  "p"' 

L=NaN; 
photo=NaN ; 
return 

end 

[ L2 , R2 ] =t ext read ( ' amO l.txt', ' %f %f ' , ' header lines ' , 1 ) ; 


[a  b] =max (photo) ; 
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R2norm=R2/max (R2) ; 
photonorm=photo/ a ; 

plot (L, photonorm, L2 , R2norm, L (b) , photonorm (b) , ' rh ' ) 

title (' Frequency  Response  vs.  AMO') 

xlabel ( 'Wavelength (urn)  '  ) 

y label  (  ' Normalized  Photogeneration ' ) 

axis ([.22,3,0,1.1]) 

legend ( ' Frequency  Response ' , ' AMO ' ) 

L(b) 


F.  TIME.M 


function  [T]=time(s) 

%This  funcion  takes  an  input  time,  's', 

%in  seconds  and  returns  time  in  the  format 
%of  hr : min : sec 

o. 

o 

% (c)  2004  B.P.  Davenport 

hr= floor (s/3600) ; 
min=floor (s/60-hr*60)  ; 
sec=floor (s-min*60-hr*3600)  ; 

if (hr<10 ) 

hr=sprintf ( ' 0%u '  ,  hr )  ; 

else 

hr=sprintf ( ' %u ' , hr) ; 

end 

if (minclO ) 

min=sprintf ( ' 0%u ' , min) ; 

else 

min=sprintf ( ' %u '  ,  min)  ; 

end 

if (sec<10) 

se c= sprint f ( ' 0%u ' , sec) ; 

else 

se c= sprint f ( ' %u ' , sec) ; 

end 

T=sprintf ( ' %s :%s : %s ' ,hr,min,  sec)  ; 
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